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The K-Shell Photoelectric Effect in the Forward and 


Backward Directions 


Arne Reitan 


d 


A 


to the incident photon. 


It is well known that to the lowest order in 
„Z, where Z is the atomic number of the element 
n question, and a=1/137, the cross section for 
»hotoelectric emission of electrons from the K- 
ell of an atom is zero in the direction of the 
cident photon as well as in the opposite direc- 
ion. This fact was first established by Sauter 
(1931). Various authors (Sauterand Wüster 
1955), Banerjee (1958), Gavrila (1959)) 
aave attempted to calculate correction terms to 
e Sauter cross section, but, at least as far as the 
forward and backward directions are concerned, 
hot very successfully. After the completion of 
the present calculations it came to the author's 
notice that Nagel (1960) had made extensive 
ocior studies of the K-shell photoelectric ef- 
ect, considering also the forward direction in 
some detail. His results, which were obtained by 
L different method from the one applied here, 
we in agreement with ours, but are not given 
| plicity for the limiting cases of emission angles 
and zt. The photoelectric cross section in the 
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Abstract. 


The first non-vanishing correction term to the Sauter cross section is calculated, by the successive 
Botn approximation method, for photoelectric emission in the directions parallel and antiparallel 


forward direction has been calculated to all orders 
in aZ for the extreme relativistic case by Mork 
and Olsen (1961). 

We are in the present paper mainly interested 
in intermediate photon energies, i. e. energies of 
the same order of magnitude as the electron rest 
energy, where relativistic effects are important, 
but where none of the simplifications that are 
possible when the photon energy is much larger 
than the electron rest energy, can be made. Since 
we are using an expansion in aZ, our results can 
only be expected to be valid to a good approxima- 
tion for small oZ. 

The matrix element of the photoelectric effect 
may, conveniently normalized, be written 


M = (Em*/m) (Zkl2)* f V3 (rase v, (n) dx, (1) 


where v, and v, are the wave functions for the 
electron in its initial and final state, m is the 
electron mass, k and € are the momentum and 
polarization vector of the incident photon, E is 
the energy of the electron in its final state, and 
0 is the usual 3-component Dirac matrix. We are 
here not interested in correction terms to the 


114 


matrix element of higher relative order than aZ 
as compared with the main terms. Accordingly 
the initial-state wave function may, in the split 
representation used by Olsen (1958), be writ- 
ten 


(2ar) (^f Eo 
WIPE AGE 6 2 
where 
a=aZm, 6 ,—grlr, (3) 


and where 7, is the Pauli spinor for the initial 
state. We have here put 71 —c— 1. 

For the final-state wave function, in order to 
obtain the results to the desired order of approxi- 
mation, we have to use the second-order Born 
approximation expression - 


V, (r) = v$ (r 3- và (r) yz (r), (4) 


where 
gą (r) — eo, 
yi()eez(—-iav-cgmd E) Î G(r—r) 
("rei ans, 
p3 — (e 2()- ia V+ Gm + E) f G(—r) 
(7 hr)(-iav'+8m-+ E)G(r —r") 
(eee fA NCE AZ SU Eu Ke, (5) 
in which 
Gr—r)=(2 ny? | (عسع لړ‎ (à — p! —żój! نل‎ s, 
1 
GA RS s 


Here p is the momentum of the emitted electron, 
B is the usual Dirac matrix, and x, is the Pauli 
spinor for the final state. In the expression for 
G(r-r’) the integration over s is as usual to be- 
performed as a complex integration, where one 
afterwards takes the limit 8->0, As in similar 
cases, we have to consider, in es course of das 


‘exp br), expecting only ph si- 
quantities to بيان‎ 
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element in the lowest order. Similary, w3 gives 
contributions of the relative order aZ to the matrix 
element, as compared with the lowest- order terms, 
and therefore has to be included in the press 


approximation. 
The matrix element (1) may in our approxi- 
mation be split into two parts 


M — M*4- MI, UE 
where 

M — (Em) (Ek/ 28 11 [LE + m) > (a p) (oe) 1 

+ z(a/2m)(o &) (o 13) 
+(E-m(E+m)1(6p(o8)1! + (o 13) (o e)] 4, , 
M'= a ZUEm[n) (Ek/ 2)? X 3 [(4/2) (E+ m (op)I; 
+ G/2) E 4- m) (o €) (oH!) 
+ (E — m? (E + my? (o p) (o e) I? 
+(2+ m) (o IŻ) (6 e) + (E+ m) (op) I$ (oe) — 
+(£— m) (612) (68)] %4. (8) | 


Here M? contains the whole contribution to the 
matrix element in the lowest order. 


The quantities I, — where the superscript j 
indicates that the term in question arises from the . 
part vZ of v, — occurring in eq. (8) are the 
following: E 


I? = (a Z71 | d? x e'2*-47(2 gr)” (a ZY*h, 
I9 = (a Z! f 43 x e/2*- 7 (r/r), 
_q=k—p, 
n —[4(1) A, B 
I: =[2(l)s 7, 
RB )68( رم‎ 
11 — [40)6/7 Z, 
dba =d? Sia! di, 


Jac d oi RUM - 
"R 
Ey 
oi: eed, 3 
Pie ود‎ 


NS 


aoa gia ag RAT A 
Fi, = Jr, Peat a, SS. 24) (amp 
(p? — ود وور در دو‎ (9) 


ORO 
It should be noted that in our approximation the 
factor (2arr 42"? 
the integral I9, 

The integrals J? and I9 may be evaluated im- 
mediately: 


(riso dtes Dries ty 


has to be considered only in 


I = 8r m (g* + a2)? [1 + (1/8) (ag| m?)] , 
19 = (i/a) 8n m (g? + a?) ? q (10) 


The remaining integrals may be Er E in terms 
pf the quantities 


I (a, b)=J a (1) 
J(a, b) =f dsr A, 


-1 7 , 
۸ د‎ „A 


where in our approximation 


11 =— (d/da) I (a, b), 

I} = — (d/da) J (a,b)=k 11 av, 1(0, b), 
I} = — i [ø J (0, 2)] (o £) (o VE), 

I} = — ży; 7(0,3), 

ROM (O) 


= — (2mk) 117 (0, b), 
I3 = — (2mk)* [o J (0, 2] (ok), 


I? = — (2mk)' J (0, b). (12) 


In the expression for I} the operator v; operates 
on J, the quantity (g vx) having been written at 
rhe end of the expression only to preserve the 
right order of the scalar products containing 6 
matrices. 

Evaluating the integrals (11) and considering, 
as far as J, is concerned, only the limiting case 
7—0, we obtain 


I (a, b) = Qn/p) (9? + à) 
| (Arctg[(£ — p? + a*)/2pa] 
=i In {(b/ 2p) (g? + ارو دوو دو — اس‎ 
J(0o,) — $p-- nk, (13) 
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where 


Re (8$) DT + [p? # - (pk)?] 
+ (2/4) (k? - pk)]), 


-1 [(p k) 


Im (8) (n/p)(247?(Q —1)— [p* £? — (pk)*]  [p&/n$ 
MERE 
Re (n) = (Ip) [p^ & — (pk! (—p 


- در 


Im (n) = (rr |p) [p? &? — (pk?) (($/4) (p k) ln ؤ‎ 
+ p* n), 
S=(p—k) (pray, 9—g7(p —P) (14) 
in which 
Q— 1- In (bm EDT). (15) 


Using eqs. (12), (13) and (14) we then ob- 
tain 
 [(mk)-* + (anl pe?) 
+ (4r ailpg?) [— (El2m? £) + (Q/4?)] , 
I: = — 27 47 [(k/mk) + (ax p/pg?)] 
+ (2r ailpg?) [ — (E + m) (k/m? k) 
+ (2p Q/9’)], 
i (o p) 14 = (p 6, — Fv — n) (ep) + to (P? S; 
+ (£& + ny) (pk) + qe +4] (e X p) 
+ (f + np) (e p) (p x k) 
+ ($+ — n,) [p (e X k)] p}, 
I; = (22? q/pg') — (2m ilp g’) [(k/mk)+ Qq Q/4*)], 
I? = (x|mpk4?) (n/2)-- 4(1 — اا‎ 
I; = (rk/mpką?) [r/2)- i(1 — QJ], 
(o p) 13 (o e) = — (2m £) V? (e p) 
+0 [£^ (p Xe) - P? £(k x e)]}, 


Ii =—2a4 


I? = —(2m E) (£p * k). (16) 
We have here written 
v E= p + £k, vin nP +k, 67) 


where 


Re (Ss) =(n*/p){— 297*+ [P € — (pk)*]-"[1 
+ (plg*) C— p? + p k)] 
+ 2 (pk) [p* £* — (عام)‎ 2 [(p k) 
+ (lg) (k? — pk)];, 
Im (Sp) = (lp) (4 Q/g^) 
+ [pg £* — (pk)]-'[— (2 pki’) 
— in ð] + 2 (pk) [p? — (pk):]-* 
[—pkln $ — (pk)/n 0]j, 
Re ($) = (n?/p) )2 ١ 
Hp? £* — (pk)?]-' (p/q?) (— £ + pk) 
T 2(p k) [p° k” - (pk) 
[P + (pla) ($* — pk)]), 
Im (Se) = (np) à — (4 Q/q*) — 4.97? (p? — £) 
+ [p E (PEC pk 
+2 (p? + pk) (p? — کم‎ 
— (ens) 22° [p k? 
— (pł?) [pk In E+ (pk) /n 6]), 
Re(n,) = Re (+) — (2/545), 
Im (n4) (n/p) || تد‎ k? — (pk)?]-* 


(Q 5/79) + (p/k) In $ 
+ 2 (p k) [7 k? — (pk)?]-* ((4/2) (pk) ما‎ ٤ 
اب رس دو‎ 

Re (ne) = (n?/2) {[p? شب‎ 


— (pk)*]-' (2/75) (— £? + pk) 
+293 [2 k? — (pk)? دم‎ 
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As expected, the terms containing the quantity 
Q disappear from the final results. 

We now introduce a coordinate system in which 
the direction of the incident photon is chosen as - 
z-axis, and in which 6 and q are the usual angular 
coordinates for the vector p. Furthermore, the co- 
ordinate system is chosen so that the polarization | 
vector of the incident photon is 


e= A(re,cFibe,) A-(r--5* 5, . (19) 


in which 
+1 for i - handed elliptic polarization, 


0 for linear polarization. (20) | 


It is convenient to write, instead of eq. (7), 
M= Mo + M, (21) 


where M, contains the contribution to M from 
Ij and the lowest-order parts of I9, I! and Ii, 
while M, contains the contributions. from those 
parts of the latter integrals which are propor- 
tional to aZ, plus the whole of M!. The lowest- 
order contribution to the cross section, summed 
over final spins and averaged over initial Es 
is then proportional to - 


«(69521 8 


19 2 


which with our Lowe: gives. 


` sin? 6 
E 0 (8, Q)= KEDI 


w :ویب‎ 
B Xx ur 


A= Z5 a9? (mlk) 8°, (24) 


rhere A. is the Compton wave length of the 
lectron. 
The first-order correction to the Sauter cross 
ction is 


1 
9, (9, P= >>) (Mo Mi+ Mi Mi), (25) 
11, x2 
or which we obtain the expression 
o, (6,p)=naZ[G, (6, p) + Gg (0, p) 

TU (8, v)IB], (26) 

„here 
1 sin? 0 a Z 1 
G, (6, Q)= : 


2 (1 — £ وی‎ 
[2 y (y — 1) (1 — B cos 6)]7* 


ES denn cos 2 9), 


— ô? 

r? + 6? 
Wj. indu ie uo Aeg 

Ga o=) (1 —$ cos 6)? 

17 ero ua لل و‎ 


y(y —1)(1— cos 9) — 1] +(y —1)(cos@—8)+ 8 
-127(7 —1)(1 —$8cos8)]-*:[—y(1—8cos6)-- 1] 


j تم‎ — 62 ۱ 
— cos | UP gs 005 2 9) |. (27) 
-utting here ۵-0, we obtain the result of G a v- 
ila (1959), which was, however, given in a 
omewhat more complicated form. This is prob- 
bly what caused that author to state that the 
irst-order correction e, does moż vanish in the 
orward and backward directions. Taking, how- 
ver cos 0— 1 in the expression for G, and 
onsidering the relation between f and y, we see 
hat و6‎ in the same way as G,, does indeed 
anish in both these cases.* 


* This fact has later also been observed by Gavrila 
(1961). 


— 
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It might be argued that it is obvious that the 
first-order correction term to the cross section 
should vanish in the forward and backward di- 
rections, because the contribution M, to the matrix 
element is zero for 0 —0, m, and so the products 
M,Mi and M,M! must also be zero for these 
angles. It should be remembered, however, that 
the expression for M, contains quantities diverging 
in the limit sin 6-0, therefore the above argu- 
ment is only valid if it can be shown that these 
quantities are cancelled in such a way that the net 
result is finite also for 60-0, m. This is indeed 
the case, and we find that the contribution M, 
to the matrix element in the limit sin 0—0 ap- 
proaches the value 
M, = r a Z (Em/p) (Eb 2)? ۹م‎ A X, (m, + ima) 

to(—zde,+re,)%,, (28) 


where 
mı = kp (E + m) [2 p+(— 2 E+ m) cos 8], (29) 


7, = (g°ln m) (E— 2m + p cos 6) [1 +(4?/2 kp) In 6]. 


Considering only the forward and backward di- 
rections, we see that the second-order correction 
to the cross section becomes 


M, P, (30)‏ وشو رل 

11, 2 
since it can be shown that for these directions 
the second-order cross term MoM} + M, Mh, where 
M, is the second-order correction term to the 
matrix element, vanishes in the same way as the 
first-order cross term (25). From eq. (30) we 
obtain 


1 
(y PE (7 — 2(1— 8 cos 7 


1 3 1 4y+1 
تلو کت س بے‎ by څې‎ tt Da 
( 27 ا‎ psp 


Jag 6 cos 9) ln [(1-- 8) 7 | (31) 


1 
0,(6,9)=q Gra Z)” 


We are here especially interested in the quan- 
tities 
(32) 


in which cpa, is the maximum differential cross 


269,50 (0)/8 max و‎ Aq FF 035 (r9 max , 
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section, the various cross sections having been 
integrated over the angle o. Taking omax as the 
value obtained when only the Sauter expression 
(23) is considered, we obtain 


pio n? y 1 
0 لس پل‎ y سه‎ (1 BE 


Kn = 


eee! 
ZAWSZE 


2 

hie (6! 63)‏ 8-87 يی 
where 0, is a solution of the equation‏ 

sin? P= 5 (cos 6, — 6) = 
D 1 
Jos yer RES z] 

iuf A PAD 
| —1 170 11-48 cs | $ SA 
and 6 
B= (3/4)? (r? + 4) (a Z}, C= (3/49 (e Z)’. (35) 


E 


we have 


1 ون بو نو‎ 7 
| -— r » 2 u re f 
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It is seen that in the high-energy limit y>>1 


the low-energy limit B<<1 we obtain 
NETA I EIL 


that the proper limits of his expressions be taken 
for sin 0—0. A comparison with the numerical 
result of Sauter and Wüster (1955) for 
uranium at the photon energy 1.33 MeV shows 
that there is no agreement between the preser 
results and theirs. However, as mentioned by 
Nagel (1960), this is due to the fact that 
Sauter and Wiister (1955) considered only 


02 


07 


Fig. 1. The ratio xo of forward to maximum intens ty 
in units of B, as a function of energy. Fully drawn li 
Theoretical curve. Dotted line: Experimental curve « 
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g. 3. The theoretical ratio xo of forward to maximum 
tensity, in units of B, as a function of energy for the 
region y = 1.0/1 و‎ 7 
art of the first-order matrix element when cal- 
ilating the cross section in the forward direction. 
Since the present theory employs an expansion 
. powers of aZ it is obvious that the results 
otained here cannot be applied to heavy elements 
ath any satisfactory degree of accuracy; the less 
| because the results indicate that the actual ex- 
insion coeffecient is maZ rather than aZ. One 
ay hope, however, that when the ratio مو "اوه‎ 
used to obtain the quantities x, the results may 
ive a somewhat larger range of applicability than 
je expressions for e, and oma, separately. 
In Figs. 1 and 2 the quantities x,/B and x4/C 
ven by eq. (33) have been plotted as functions 
y. The curves qualitatively show the same be- 


و 


haviour as that found experimentally by Hult- 
berg (1959) for Z=92, but, as expected, the 
quantitative agreement is not very good. The ex- 
perimental points obtained by that author for x, 
are plotted in the diagram of Fig. 1. It is seen 
that there is a fair agreement between the theore- 
tical and experimental results for x, in the region 
around y=2. For higher energies the theoretical 
curve lies above the experimental one, whereas 
the opposite is the case for lower energies. 

In the case of x4 the experimental points ob- 
tained by Hultberg (1959) for Z—92 have 
not been marked in the diagram, since the theore- 
tical results are here 10—100 times smaller than 
the experimental ones. 

In Fig. 3 the quantity x, has been plotted 
as a function of y for the low-energy region 
y=1.0—1.3. It is seen that x,, as given by eq. 
(33), vanishes, not only for y—1, but also for 
y=5/4. 
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Abstract 


In the present paper we shall be concerned with an infinite bilinear series of Gegenbauer poly- 
nomials. In the first section a sum formula will be found, and the convergence will be discussed. 
In the next section a special case of the sum formula will be applied to Schrédinger’s energy equa- 
tion in momentum space for a particle in a Coulomb field. The appendix contains a lemma on the 


convergence of a power series on the circle of convergence. 


I. Sum formula 


We consider the infinite sum: 


c= —1 .v v 
S=ZEm![T(n+- 27(( Cele) GZ, (1) 
x and y are real numbers in the closed interval 
[—1,1]. v is real and positive. z is complex, |z] 
rk 
The symbols C,(x), 7» = 0, 1, 2, .. ., denote a 
set of Gegenbauer polynomials. These can be de- 
fined by means of the expansion 


OO 
(1— 25x ړت‎ 2 G„(2)9*, (2) 


whose radius of convergence is |b|=1, when 


—1 <x <1. 
In order to examine the convergence of the 


| series (1), let $y denote the N'th partial sum. 


| 
Í 


| 


Then we obviously have: 


22 


2 
max 


| Sl ZE [P (nb 27]! i (x) 


where | C; (X) |max is the maximum of |C} (x)| in 
the interval [—1, 4-1]. According to Bateman 
(1953), vol. II, p. 206, one has 


v | د‎ T + 2v) 
| (x) max = Cx (1)= I (2v)n! 


Hence: 


[uut pomeS a ve pe 
| 1 
|< کې‎ Šį! T Jeb. 


Since we have assumed |z| <1, this demonstrates 
that the series $ is convergent. 

Our first step in summing the series is to make 
use of the formula 


z! [T (v) G, (x) Cy) = 
1 
T (n + 2v) 27% Je (r)1-2)7!d, (3) 
-1 


z xy—t(1— × (1 — y^, 
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which can be found in trigonometric notation in 
Bateman (1953), vol. I, p. 177. Then we get 


= 1 


$— 2v (r (v? Vy ZI (1—By-1dr. (4) 


21—0 


-1 


Here we can interchange the order of summation 
and integration. This is true because the series 


oo 


3 e) |z|<1, 
n=o 

is uniformly convergent in the domain of integra- 
tion. (One should notice that when ! is in the 
interval [— 1, +1], so 1s also z. Thus in the 
domain of integration: c e ce Then 
the uniformity is easily seen by means of the Weier- 
strass M-test.) 


Hence: 
1 oo 
Sc OS at Ah Dy Gal testas 43) 


-1 


The last sum is the generating function, 


1—2zv4 zy =$ 6 (z) 2”, 


+ = o 


Consequently, 
See Lv A pee 
1 ۶ i 
jaa — fy) (1 DT (6) 
E 


We introduce the complex quantity, 
یزور کا‎ z) 1—)(1— O) 

which, for |z| < 1, takes all values in ihe complex R 
plane, d eer on the line joining — 1 and 
--1. Then the generating function can | 


- ten as nnm 
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— According to Bateman (1953) vol. I, p. 155, one 


where (= } (z+ z! — 2x) (1 — 3) (1 قر‎ 


hand side, we put — 7 < arg z € 7. Then the 

functions z” and ($+) will have discontinu- 

ities when z passes the negative real axis. Their 

product, however, is continous (when |z| < 1). 
In this notation we get the following sum: 


څې هوس د rpi- y‏ )ماو سو 
1 
oe Jo — AYA n dr (9)‏ 
1- 


has the following integral representation of Legen- 
dre functions of the second kind: 


1 
QD zaj 1-2) (GED dt. (10) 
=! 


This is valid when Re v> 1, |arg (£7 1)| <<a 
and ¢ is not on the real axis between — 1 and +-1. 
Thus we have the result: 


Sa Dy mF (n+ 2v) C G9 Ch (y) 2" 


= 2-2 [ro (1 -x 20) 2. | 


2 Qa notans AJ 


A 


and |z | 1. =" 


- Q„(6) is a single-valued and regular function i 
the ¢-plane when cut along the ai a (= 


مه + 
کا( 


| 


— ©. It has logarithmic singularities in (= 
and also a singularity in ¢=oo when v ېې‎ 
Since S is a power series in z, the expt 

4, ( is an analyti Funct: Wł: 


EW 


S= [I (2 y) (12 2% 4-24) 
Fata لاد‎ 


| If we apply the trigonometric notation, 


(13) 


X—cos a , y=cos$, 0< a D « m, 


and introduce the orthonormalized set of func- 
tions, 


c, (a) = 
(n+ v) n! 


E رت‎ 


2 
٢ ext a)d a =0,,, (15) 
0 
q. (11) can be written briefly as ^ 
"v 7 
UN c, (a) cz (B) eant اا‎ )16( 
و لس‎ s 
1 m 
e ger ca ae ‘os 6 
„ao E oi cs (17) 
sin a sin Û 
We also have the interesting form: 
zr k 
Qa) n Tes — cn (a) cx (6) (18) 
Let us now consider the case 7 |=1. Our starting 


oint is the function 
s(zj=z"! سه‎ 


hich is the sum of the series (16) when |z] 

1. In this domain s(z) is an analytic function 
f z. sina and sin 8 will be assumed to be different 
rom zero. Otherwise all terms of the power series 
re identically vanishing. s(z) has generally four 
ints of singularity when |z|= 1, namely when 
Ee 


z= e, cos P = cos (a F p). 


Let z, be one of these points and let z be a point 
in the neighbourhood of م‎ Then: 


jug E ا‎ 
i JA oom) 
2sinesinf 7 > 
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As can be seen from the integral representation 
(10), the singularities of O (€) in € = 1 and 
S=— 1 are of logarithmic nature. The divergent 
part is proportional to log (¢—1) and to 
log (¢+1), respectively. Thus, in the neighbour- 
hood of the singular point z, on the unite circle 
the divergent part of s(z) is proportional to 
log (z — zo), 


i.e. |s(z) ~ —log|z—2,|. 


In the neighbourhood of the singular point z,— 
e's, we can always assign a positive integrabel 


function m(q). say m(q) = |P — Po : such that 
s(pe?r) < m (p) far pai Consequently, it is 
possible to assign a positive integrable function 
M (p) such that for any ۷ and any p< 1, we 
have |s(z)| < M (9). 

From eq. (17) we get the limit, 


COS P — cos وو‎ D (19) 


lim s(z)==e"*9 Qy.4 
sin a sin f 


p—>1- 
The symbol, = 7:0, is used to define the correct 
function in the case that the argument jn OE 
on the cut from — oo to 1. The upper or lower 
sign is valid according as م‎ 0. When q is equal 
to zero, the argument is larger than 1 ifa-F$, 
and equal to 1 if a =8. When Q is equal to 7r, the 
limit can be found by noting that s(z) is invariant 
under the transformation P >P — 7t, «— 7t — a. In 
either point, p=0, rr, the limit is a continous func- 
tion of ~ (except when the point represents a sin- 
gularity). 

In any interval of p, which does not contain any 
of the singular points, the limit (19) is obviously 
a continous function of bounded variation. 

Refering to the appendix, we can conclude that 
the power series (16) is convergent on the unit 
circle when cos p -- cos(« = 6). Ihes sum is given 


by 


oo 


2) ae c, (a) c; (B) e"? 


=e? Qva (69 4 2:0), (20) 
__ cos — cos a cos B 
where $= BAY. , (21) 


upper or lower sign according as p20. 
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For p=0 and P= r we have, 


ve TL v v ak. 
> nd-v Cn (a) Cn (6) 22 e| 


1-0 


= os 
1 Ed (22) 
sin a sin Û 


and 
= 7 
چو کا‎ cn (a) cn (P) 
A | 1 + cos a Cos f د‎ 
em سه‎ sin a sin f } (23) 
respectively. 


Writting eq. (20) in the form, 


oo 
TU v v : 
i (n+) p 
den c, (a) Cn (Be 


= Opi (Gy —2 0) P>; 


and combining this equation with the complex 
conjugate one, we get the results: 


> 7 
Duy 


n=o 


y i [o (QI Tees MEN o] 


c, (a) c, (B) cos(a + ») p 


= Quy (ea (2 4) 


2) > e (0) c; (8) sin (n + v) p 


— mine 
- [06 )ولب )0 .1۰ ده‎ 


M 
DE TZ 


= Pa (o) 


Erik Eriksen 


for a particle with mass m and energy E <0 in 
a Coulomb field V(r) = Ze?/r. Unites will be 
chosen such that h =1, i.e. p*=— NE 

In momentum space the equation is: 


(& E)200+ 


2m 
-3 
(2 mr) ‘| U(|k' —k|) 2(k)2k' =o. )27( 
Q (k) and U(k) are the Fourier transforms: 


0 (عل)‎ am) |r cok dr, (28) 


pe an) Vine dr = 


35 


— (27) * Am Ze kt. 


The energy-eigenfunctions in momentum space 
have been found by various authours applying 
different methods. Pauling and Podolsky (1929) 
made use of a Fourier transform of the well known 
solutions in coordinate space. Later Hylleraas 
(1932) deduced a second order differential equa- 
tion i momentum space, which he was able to 
solve. Fock (1935) solved the integral equation 
in momentum space, eq. (27) above, by a trans- 
formation of the k-space to a hypersphere in a 


by spherical harmonics. — 
bx hy shes? Tni vedi Fate 
SETS Plecy sits tis emia Me 97 


node 


e substitute the expression (31) for @ (k) in 
he integral equation and carry out the integration 
ver the angles. The result is the following inte- 
gral equation for the radial part: 


E n" z) MO fe (4,5) frlds=0, (33) 


2m : 
0 


here the kernel — in the general case of a 


central field V (7) — is given by, 


K,(k,)= I rV (r) (sk)? Joa Gr) Jip (kr) dr. (34) 

or the Coulomb field: 

z (k, 5) = — ze CWE (sr) J:+4-(kr) dr. (35) 
0 


he integral is convergent except when > 
According to Bateman (1954) voll, p. 50, 
the kernel is: E 


2 2 2 
Kit, j=— 2 o, m 


2h (56) 


ze | gauge 
ته یئ‎ o (E 7 o di. (37) 
0 


his integral equation can be converted to a 
ymmetric one in different ways. Introducing the 


rguments œ and f, defined by 


a 
—(— 2m E)? ee s=(— 2m E)” «^ (38) 


(39) 


'we get the equation 


1 — cos a cos B 
sin & sin B 


v2 | * o( ٨ 48, )40( 
| 0 
where A= Ze? ) 21 (41) 


This equation can be solved by means of the sum 
formula (22), which resolves the kernel in a 
bilinear series of orthonormalized functions: 
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1/25 
! o 1 — cos a cos f | m 
71 sin a sin Û 
A et 
7-81 73-1 
z - - 8 A Ą 
FTES (62 0. (42) 
Consequently, eq. (40) has the solutions, 
gi (e) c T (a),n=0,1,2,:*-. (43) 
with the corresponding eigenvalues, 
A= Ze ( —m/2E)^ — n-- 14 1. (44) 


Appendix. 


We shall prove the following lemma: 


OO 
Let f(z) — Zi c, 2” be a power series of unite 
radius of convergence, and let there exist a posi- 
tive integrable function, M(q), such that 


| f (pe?) | < MQ) 


for e <1 and every q. 
Let us suppose the Abel sum, 


ES : 
Flp) = lim So, p” e", 
0—1— 2—o0 
to exist everywhere or almost everywhere on the 
unite circle. Let a be an interior point of any inter- 
val (@,, p) within which F(q) is of bounded 
variation over the set of points where F(q) exists. 


oo 
Then, Z(«) and the sum 2c, e^ exist and 
1-0 


ate equal. 


Proof. 
|F(g)| is defined everywhere or almost every- 
where on the unite circle. Due to the existence of 
the function M(g), |F(g)| is integrabel, and we 
have: 


Tat) dy fira 


—u 


7 
lim Î | f (pe?) dp. 
(See E. W. Hobson: «The Theory of Functions 
of a Real Variable», vol. II, p. 323, Dover Publi- 
cations, New York, 1957.) 
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Consider the Fourier series associated with F(q), 


FA) ~} a, + (an cos np + b, sin n9), 


کد 


IU m 
ır a, =) F(g)cosng dp, 7 b, = f F(p)sin ng do. 
—T پټ‎ 

The Fourier coefficients are independent of the 
values we assign to F(q) in the set of points 
where F(q) is non-existing. This is true because 
these points are a set of measure zero. Let a be an 
interior point of any interval (pı, 9») within 
which F(p) is of bounded variation over the set 
of points where F(q) exists. Then, according to 
Fourier's theorem the series (45) is convergent 
for p=a. 

Now we shall show that the Fourier series (45) 
actually is the power series f(z) = Zc, z" for z=. 

Due to the existence of the function M(q), we 
have: 


m F(q)cos np dq = J c: f (pe) cos np dp = 


=lim i f (pe^) cos np dp = 


0—-1— —x 


= lim n Z c, p" CET cos np dq. 

0—>1- —a n=o 
When e < 1, the series in the last integral is uni- 
formly convergent, and the integration can he. 
carried out term by term. The result is: 


(* د‎ 
Atte JH 


„kom 
4874] 


E 1 
وم‎ uM Smile په‎ 
+ Ba AS j 


r» boaMAU "ې‎ 
at | hi Hn ndt "an 
[a SL. 


OC atsi " "E خم‎ 


wo es 


va Las Fani له‎ EXE 


کې 
u^ :‏ 
.ما 
46871 


د پخ 


Erik Eriksen 


= à "x 
BE 


Í 2Co, 0 
a, = 
‘i | (ago GIŻE 
Quite analogously we find for the Fourier coeffi- 
cient by, 


b, Zi, N=1,2,°° 
Thus, the Fourier series associated with Z (p) is 


oo 
Z c, e”? , and it is convergent for p = a. 


n=o 


According to Abel's theorem on continuity up 1 
to the circle of convergence, we have 


oo O oo H 
Ze, e"* —]im De, p" e/"* — F(a), 
ور‎ 0—1- n—o RN 


Hence, F(a) and the sum 2c, e" exist and 


n=o 


are equal. 
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g¢ ^ Abstract 
P 

The derivation of the formulae for neutron scattering by homonuclear diatomic molecules using 
the Fermi pseudopotential in conjunction with the Born approximation is outlined. The quantum 
mechanical state of the system before and after the interaction is taken into account in a detailed 
manner. 

Numerical results are presented for Hy, Dı, Ng and O, and compared with available experi- 
mental data. The properties of H, and D, as cold moderators are discussed briefly. 

Two possible sets of values for the spin dependent scattering lengths of nitrogen are obtained 
by fitting the theoretical curves to the experimental data, 


2a, +a_= (2472-0023) 10-8 cm 
a, — کے‎ t (0.73 £0.10) 10-12 cm 


The magnetic scattering from the Os-molecule is determined by subtracting the calculated nuc- 


1. Introduction. 


The slow neutron wavelength is of the same 
rder of magnitude as the internuclear distance of 
diatomic molecule. One therefore gets inter- 
ference between the waves scattered from the two 
uclei. Spin effects can give different amplitudes 
to the two waves, in some cases also different 
hases. 

If the rotational or the vibrational motion of the 
molecule is changed during the scattering process, 
'we call the scattering inelastic. If they remain the 
'same the scattering is called quasi elastic, only 
translational energy can then be exchanged. 

Use of the Born approximation in conjunction 


— * Work sponsored jointly by Institutt for Atomenergi, 
Norway, and Reactor Centrum Nederland, the Nether- 
fands. 


lear scattering from the experimental values of the total cross section. 


with the Fermi pseudopotential makes it possible 
to find the transition probability between two 
quantum states. By averaging over all initial states 
and summing over the final ones, the total scatter- 
ing cross section can be obtained. | 

This method was developed by Schwinger 
and Teller (1937) (later referred to as ST) 
for Hs and by Hamermesh and Schwin- 
ger (1946) (later referred to as HS) for Do. 
Their main purpose was to determine the spin 
dependence of the neutron-nucleon interaction and 
numerical calculations were only performed for 
very low neutron energies and for the lowest 
molecular states. 

Brimberg (1956) has calculated the cross- 
section for Hg at room temperature by this method 
and also corrected for zero point vibrations. 


"^ real هه‎ ii base 
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Later much work has been done to find simpler 
and semiclassical approximations. With the com- 
puting machines available to-day, the case of 
diatomic molecules can easily be handled by the 
procedure indicated above, which should give the 
most detailed and accurate information. 

In Section 2 we develop the formulas for neu- 
tron scattering from a homonuclear diatomic mole- 
cule in general and discuss various approximations 
for the molecular wave functions focussing our 
attention on the four molecules Ho, Dg, No and 
Os. The vibrational effects on the scattering are 
considered. They are found to be of little impor- 
tance for the two heavier molecules, but should 
be taken into account for Hy and Dy. At room 
temperature these molecules will be in a thermal 
distribution among the lowest rotational states, 
excitation of vibrations being unimportant. 

In Section 3 the procedure for numerical evalu- 
ation is outlined. In Section 4 graphs of the neu- 
tron scattering cross-section for Hz, Dy, Na and 
O, are presented and compared with experimental 
data, covering the energy range up to the thres- 
hold for excitation of the second vibrational level. 

In the last section the limitations of the met- 
hods applied are briefly discussed. 


2. Scattering formulas. 


The differential cross-section for neutron scat- 
tering by a molecule of mass 2 Am will in 6 
Born approximation be given by 2 


b REAR" 
i esz 2«53] 7^ 


x| <w, Uv; |20 (1) 
Subscripts and f indicate initial and final states 
respectively. p is neutron momentum, mm neutron 
mass, V normalization volume, ¥ the wave func- 
tion for the whole system, and dQ an element of 

solid angle. The interaction operator U will in 
general be ap a دل‎ For a single nucleus of 
mass Am and spin s, it can be written 
dia DISE MAS مهه‎ uli to sombnogeb 

wi نه‎ Hah. x جار‎ lata 
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as given by HS. Here 8 and o/2 are the spin - 
operators of the nucleus and the neutron, a, and 
a. the scattering lengths for parallel and anti- 
parallel spins respectively. 8 (ra — r) is the Dirac 
delta function for the position vectors. s 
The total scattering cross-section for a free 


nucleus will be given by 
(3). 
If the scattering is spin independent this reduces 
to a=4ra2. In the case of scattering of point 


particles on hard spheres, one sees that the scat- 
tering length will be equal to the radius of the 


sphere. 


2.1 Cross-sections for a homonuclear diatomic 
molecule. 


The total interaction operator is then most con- 
veniently divided into the sum of a symmetric and 
an antisymmetric part, 


eo ale G+ 1) ay + 25a 
+ (a, —a.) o S] [8 (m — rı) + 8 (ra — r2)]} (4) 
gue — pes po 
x[Bm—n)-85(—r) (5) 
where i 
l S=sits, 6ه‎ 


Y ; 292 
| EU 
1 be of the form ZA : 


tric part g 


_ The we 


71 the rotational quantum numbers. r = Fı — ro 
ji the internuclear distance. 

| Since we are treating homonuclear molecules, 
jruantum statistics impose the following relation 


ESET 4 (8) 


yhere e is dependent upon the electronic symmetry 
d has the value +1 for hydrogen, deuterium and 
utrogen, and — 1 for oxygen. 

We assume the energy change in the scattering 
be independent of the spin and projection quan- 
um numbers. Eq. (1) can therefore be summed 
ver the corresponding final states and averaged 
tver the initial ones. By substituting eqs. (4), (5) 
nd (7) into (1) one then obtains, by the met- 
10d of ST, for a transition from an initial state 
47 to a final state vj 


M 
+ 


a 

bf | (A+ | 2 
BOs = سي‎ = z! > 
Tue 2 1 B; EE <A 


CZ | RZE ېږ‎ EXP zq r) NZ d (r) D dQ (9) 


CZA 


here 


q= (pi — pı)/ 2h (10) 


s proportional to the momentum transfer. The 
xplicit dependence on spin quantum numbers is 
liminated because of eq. (8) but the spin depen- 
ence is contained in /he quantity B, given by 


5 


+ S(S+1)(a, — a / 5 (254-1) (11) 


»— و‎ odd: 
Bones {(20-+ ۱)2 1) a )به‎ 1) — 
5 


— 5(S+1)](a,—a_)* Wer (25+1) (12) 


| 
The prime on the summation symbol indicates 
at the values of S are restricted by eq. (8). 
Aq. (11) will of course correspond to coherent 
nd eq. (12) to incoherent scattering. 

For the cases of interest here, B,» is given 
in table I (assuming that e=1). 
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Table I 
EAD By 
O | e | e | 4ra’ 
| 
| 6 e 


wj 


|o e A == zaa 
o|o " [Be +a- + 2 (a, — a_)?] 
ele 7 [4 (22, +a_)?-+5(a, — a] 
e padla, a 
1 
o| e = (a, —a_} 
o |o 5 4 (244 ٨٢7٢ TA 


To reduce the matrix element in eq. (9) we 
will assume that the vibrational and the rotational 
part of the wave function can be separated, or 


Vo jm (£) =*,(r) Yin (9, o) (13) 


where Yj», the wave function of a rigid rotator, 
is the normalized spherical harmonics of order j, 
and ®,(7) is the radial wave function. 

Using the well-known expansion 


exp (qr) = EH V 4r (2k+ 1) j (gr) Yeo (9, p) 
where j, is the spherical Bessel function of order 
k, the angular part of the matrix element in eq. 
(9) can be expressed as a Clebsch-Gordan coef- 
ficient, and the summation over m and m’ can 
be performed utilizing the orthogonal properties 
of these coefficients. This is an alternative to the 
method used by ST. The result is 
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a v! 


A n 
يو‎ a ri 
I; 


2(2k+ 1C [J ja (qne, dur dr) dQ (Q4) 


where, by setting L= } (j ٣( +4), 
(2/ — 25)!(2/— 25')!(2/— 2%)! (7!)? 


Cy, = its; 
v = IFO ۹۱ 
if È 
de aed A W 
otherwise 
Cj 30 


We will here be interested only in the total 
scattering cross-section. The only quantity depen- 
dent on the scattering angle is g. From eq. (10) it 
follows that 

1 
28, bę— وی‎ 0]? 


)16( کو ول 


dQ — 27 sin 0d0 = 2r (45?/p; pp) a dg (17) 


In terms of the neutron energy, E, in the labor- 
atory system, we have 


2 (2mE)? (18) 
A 
ورا مر‎ s =" 2 o yi (19) 


where E,; is the rotational and vibrational energy 
of the molecules. Eq. (14) can now be inte- 
grated over the scattering angle to give the total: 
cross-section: l 
Toj uj = 8m (b/p)? X 

Q5)" (6+ py) 
>= n 1 


[total q dq 
7 (257 |p: gy] a i 
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2.2 The molecular wave functions. 

The various approximations for 91 1 
wave functions are discussed in the book by 
Herzberg (1950), Spectra of Diatomic Mole- 
cules. We shall here assume that the vibrational 
motion is harmonic. The two lowest radial wave 
functions are then 


4 


ay (0 V (aln) r exp (=3 a (r — r") (23) 


4 
V (alr) ^V 2a (r — ra) >< 


$(r)= 
Xexp( - } a (r - r?) (24). 
Here r, is the internuclear equilibrium distance and 
a= Ambv | 2b (25) 


where bv is the vibrational quantum. The rotatio-- 
nal wave functions was chosen above to be that 
of the rigid rotator. The molecular encres is 
then given n 


EID S )B — Q6) 


where B is the rotational energy constant. For a 
non vibrating molecule it is given by 


B=B,= b?/Am re rz (27) 


For a vibrating molecule the factor r shou d 
be replaced by the expectation value of 72 
(Herzberg (1950) p. 106). The numetical 
value of B used in our calculations is corrected for 
the vibrational ground state. The value of the 
vibrational quantum was chosen to be hy=. 
Ej). The most important molecular constants 


* 
Vr T 


listed in table. ال‎ Te | 
2 : ` pa m . 
hb سه‎ Ren dup prote ۶ Fe beja 
a از‎ R e M 1 
Numerical values of molecular 


هر سي 


A>>(a) ? (28) 


expressed in terms of the neutron energy E 
1« « 2« (Abv]2E) (29) 


ne Gaussian factor in the vibrational wave func- 
on will to a good approximation act as a Dirac 
Ita function, i.e. 


! (8 (r — rji 


1(8 (r— rj 
q. (21) can then be simplified to 


[E 


$, X (r— re) (rr) 


Ik = 9$. Sou "Jk (4 r.) 


٠۰3 Partial cross-sections. 
Introducing, for the ease of ENS the 


iimensionless energy variable oo 


Æ 


¢= (E/B (1 + 1/2A))? (33) 
nd the constant 
1=(rl25) ):8/)1 + 1/2A)f (34) 
e get from eq. (20) and (32) 
7 
oj =(A L)? Ay. Je (ux) dx )35( 


vhere the integration limits are given by 


Ie =r(O+(62?+7G41)—/' (7 +1)? | G6 


The subscripts referring to vibration are dropped, 
ind we will call this the rigid approximation. 
Taking the vibrational motion into account, one gets 


da 


(ALY? ZZ [i x Iova (x)] dx (37)‏ با 
the limits now M‏ 
[eid i6 0-7 9 3-0‏ 
+w - v) bv/B | (38)‏ 


| 

if is seen from condition (28) and the numerical 
alues in table II that the influence of the vibra- 
ional motion will be most important for the ligh- 
er molecules H, and Dg. Since j,(x) is an oscil- 
ating function, it can be seen from eq. (21) that 
he predominant effect of the zero point vibrations 
ill be to reduce the scattering cross-section. In 
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physical terms, the constructive interference be- 
tween the two scattered waves is partially destroyed. 
For A > (a)-% this effect will always be more 
important than the contribution from excitation of 
the higher vibrational levels because of the ortho- 
gonality of the wave functions in eq. (21) if vv. 

The total scattering cross-section for the initial 
state v, / is obtained by summing over all the final 
ones, 


(39) 


Om = Pw, uj” 
uj" 
2.4 The temperature dependence of the scattering. 
For a gas in equilibrium at a temperature T, 
the probability of finding a molecule in the state 
v,j is given by the Boltzmann distribution 


P, )1( NQjtiz ‘(25+ 1) exp( — E, £T) (40) 


The prime on the summation symbol indicates 
that the summation over $ is restricted as discus- 
sed above, The normalization constant N is defi- 
ned by: > p, 1 

v7 


The average cross-section 1s then 


o = È P; (T) 9 7j (41) 
v7 

If the thermal translational velocity z of the mole- 

cules is comparable to or larger than the neutron 

velocity v, then the effective scattering cross sec- 

tion e(v) which can be observed is determined by 

the average 


o(v)=7" f|v—ulo(|v—al)N(ujd(u) (42) 


The velocity distribution N(u) is assumed to be 
given by the Maxwellian expression 


N (u) = (Am / s ET) exp ( — Amu? / kT) (43) 
By introducing 
Am 
S = y 4T v (44) 
Am 
t y T. (v — u) (45) 
eq. (40) can be written 
o=2n 3 7 exp( — 5)» 
co 
x feo ) exp ( — £) sinh (257) dt (46) 
0 
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Dy series expansion it is seen that when s+ O 
(47) 


regardless of the energy dependence of o. The 
total cross-section which is most easily observed 
is then 


— 
o > (const.) E * 


(48) 


Ga‏ +6 چو نا 
The absorption cross-section o, will here be assu-‏ 
med to be 1/v-dependent.‏ 


3. Numerical methods. 


The numerical evaluation of the formulas deri- 
ved above has been made using the Ferranti Mer- 
cuty Computer FREDERIC at the Norwegian De- 
fense Research Establishment. 

We have computed the total scattering cross 
section for neutron energies up to 2 bv for the four 
molecules considered. For all of them the cross- 
section has been evaluated in the rigid approxi- 
mation given by eq. (35). The computer pro- 
gramme was made quite general, the only input 
needed for each molecule being the constants de- 
pending upon mass, spin, scattering length and 
temperature. 

For Hs and D, we have also calculated the 
cross-section using the oscillator model described 
above. The computation of the double integral 
involved is a very time-consuming process even 
on an electronic computer, we have therefore 
expanded the spherical Bessel function around 
rer, For the case v=0'=0 we get from eq. (21) 
and (23) by expansion to 4th order and inte- 
grating 


Los = (al)? [ja (qr) exp (= a (r— ra? ) dr 
= jg 60) + (63/48) © ( a 
+ (4/32 8%) 7 (x (49) 
where 

X= a Te 

T 
The lower integration limit has been extended to 
— w which introduces negligible errors compared 


to the other approximations involved. By using the 
recursion formula 


EE 1 ye TEŻ : 
JE" RIP iSi 


j4 (X) = ja — 


eq. (49) can be expressed as 
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look = fk t و‎ Mes 22 1) -x|+ 5x 
x [ (2 + 1) (e+ 2) (A+ 3) (6 + 4) — 


- 2(B 324 او‎ 


i x MD 2 a 
ک‎ gg ey al (50) 


For excitations into the first vibrational state we 
get in the same way from eq. (21) and (24) 


PE 
Na =a y? f (r — re) jk (Gr) expC— a (r — r2) dr = 
0 


REP | x 
- او‎ A04 


a (5) : 


2 1 3 | 
= y2\- (2 + dl 3)+ 1 


4 xt څې اد‎ E {+ 1+24+32+ 05€ 
x (k + 5)+ 2 (k? + 104 + 19k+- 30) x? — 


2 x 
—(k+ 5» بس‎ ys ig 
<|e CEG KS sA 2523 597 + 58k+- 


+120 — )2 + 2k + 20) 86) lia (51) 


4. Results. 


Graphs of the total cross section at room tempe- 
rature vs. neutron energy are in the following pre- 
sented for all four molecules and compared with 
experimental data available. We also give o; for 
j-values of importance at room temperature, which 
makes it possible to estimate the cross section at a 
lower temperature without starting from the be- 
ginning. For a few cases oyy’ is also given. a,” 
gives information about the probability for energy 
exchange in the scattering, and is therefore essen- 
tial if one wants: a detailed picture of the slowing 
down process. 

All numerical cross-section values are given in 
barns per atom. 


v 


The following values of the scattering lengths 
ete used 


84 =5.28X 10" cm 
a_ = — 23.69 X 10^? cm 


[The amplitudes for the various transitions as given 
y table I for s= 4 then becomes: B,, = 0.44 
arns, Be = 19.92 barns, Boe = 6.63 barns, Boo = 

13.73 barns. The relative smallness of B,, gives a 

characteristic feature to the curves for 9;; which 

re shown on fig. 1 as calculated in the rigid appro- 
imation for the two lowest initial states. 
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The cutves of fig. 1 are of particular interest 
since liquid hydrogen is now being used as a 
moderator to obtain very slow neutrons. We may 
assume the scattering properties to have the same 
gross features in the condensed as in the gaseous 
state. (This is indicated by Raman spectra, see e.g. 
Van Krankendonk (1950) and Bhat- 
nagaretal. (1962).) Hydrogen is in the ground- 
state almost transparent for neutrons with velocity 
below 2 km/s (E= 0,02 eV), as can be seen 
from fig. 1. It would indeed be advantageous to 
have the next lowest state present, which will 
cause further slowing-down by  quasi-elastic 
scattering, the cross section for energy absorption 
by the neutron is seen to be small in comparison. 
Now, as is well known, para hydrogen (j even) 
and ortho hydrogen (j odd) form two metastable 
modifications. If hydrogen is liquified by cooling 
from room temperature, approximately 1/4 of the 
molecules will be in the state j=0 and the rest 
in the state j=1. By adding a suitable catalyst 
this ratio can be changed, for our case a value 
around 1:1 would probably be a good choice. 
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Fig. 2. The scattering cross section for the 4 lowest 
rotational states in hydrogen gas. Broken curve: Effects 
of zero point vibrations neglected. Contribution due to 
excitation of the first vibrational state shown separately. 


Fig. 1. The cross section for a rotational transition 
from an initial state j (—0 or 1) to a final state j for 
hydrogen gas. 
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Graphs of o; are shown in fig. 2 both for the 
rigid and the oscillator model. The contribution 
from excitation of the first vibrational level is 
shown separately and not included in the main 
component. 

In fig. 3 the total average cross-section 1s 
shown and compared with experimental data by 
Melkonian (1949). An absorption cross- 
section of 0.33 barns at v—2250 m/s has been 
included, assuming 1/v-dependence. 
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Fig. 3. The total cross-section for hydrogen gas. 


4.2 Deuterium. 


From differential cross-section measurements, 
Hurstand Alcock (1951) have found 


2a, +a-= 1.346 X 107? cm 
a, — a- = 0.567 X 107? cm 


The curves for ø, are given in fig. 4 for the 
two lowest initial states. The cross-section for in- 
elastic scattering is seen to be small in the sub- 
thermal region. If liquid deuterium is used as a 
moderator, the most important slowing-down 
process will therefore be exchange of transla- 
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Fig. 4. The cross-section for a rotational transition 
from an initial state j (=0 or 1) to a final state 7’ for 
deuterium gas. 


tional energy. Since the deuterium molecule is 4 
times as heavy as the neutron, a rather large mode- 
rator volume is needed compared to what is 
necessary if light hydrogen is used. The advantage 
of deuterium is of course that the absorption cross- - 
section is negligible. 

Curves for o, as calculated by the rigid ap- 
proximation are shown in fig. 5. The correction 
for the zero point vibrations will to a good — 
approximation be the same for all initial states, 
only e; is therefore calculated using the oscillator 
model, j=2 being the most probable initial state 
at room temperature. The part of the cross section 
due to excitation of vibration is shown separately. 

A graph of the total cross-section at room - 
temperature is shown in fig. 6. All corrections - 
mentioned have been applied. 

As far as we know, no experimental data are 
available for the total neutron cross-section of — 
deuterium gas. 
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jig. 5. The scattering cross-section for the 7 lowest 

lotational states in deuterium gas. The broken curve for 

2 is corrected for zero point vibrations. Contribution 

e to excitation of the first vibrational state shown 
separately. 


43 Nitrogen. 

99.63 % of ordinary nitrogen consists of the 
sotope N7, in the present calculation effects due 
D other isotopes are neglected. 

"The spin of the N7-nucleus is 1. Nitrogen is 
erefore expected to give some incoherent scatter- 
ng, but from the data available it was not pos- 
ble to give a quantitative estimate. 

| As a first trial, we therefore assumed the scat- 
fering to be coherent and determined the ampli- 
ude from the value of 9.96 barns for the free 
itom scattering cross-section given by Melkonian. 
| In fig. 7 we show o; for j=10 which is one 
if the most probable initial states at room tempera- 
lare. For comparison we have also included the 
nost important incoherent transition at low ener- 
Lies, jï =9, giving it the same amplitude as the 
herent ones. 
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Fig. 6. The total scattering cross-section for 
deuterium gas. 


Curves for c; ate shown in fig. 8 for j=0, 
10 and 20. As the diagram indicates, e; varies 
little with j for neutron velocities above 2.5 km/s. 
In this region e,, may thus represent the room 
temperature average over all initial states to a 
sufficient accuracy. 

When the total cross-section was calculated and 
compared to experimental data, a significant dis- 
agreement was found. The incoherent cross-section 
was therefore also calculated. In fig. 9 the ave- 
raged values of the coherent and incoherent cross- 
sections are shown. To simplify the treatment, 
the same weight was given to the two coherent 
transitions (odd-odd and even-even) and to the 
two incoherent ones (odd-even and even-odd). 
The error introduced by this is negligible because 
we are averaging over 20 initial states. 

Since the energy dependence in the subthermal 
region is so different for the two curves in fig. 9, 
the magnitude of | a, — a. | can be determined 
total cross section measurements of nitrogen gas 
at room temperature. 

The experimental data available so far are not 
too accurate, but the data of Melkonian 
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Fig. 7. The cross-section for a coherent rotational 
transition from the initial state /— 10 to a final state j’ 
for nitrogen gas. j —9 included for comparison 
(see text). 


(1949) were by the least square method fitted to 
a sum of the two scattering terms shown in fig. 9 
and a 1/v-dependent absorption term. 

One then obtained 


24. lg =2 29 T 0.08) RIO em 
= د092‎ O 15) 210 76m 
a22 + 0.2) barns ai 1 2777 


The sign of the first term is known from neutron 
diffraction work. The value of the absorption cross- 
section is too high. Actually, the absorption term 
is a bit too similar in shape to a combination of 
the two other terms to make a determination by 
this method very accurate. A second fit was there- 
fore made where the absorption was subtracted 
from the experimental data. A tabulated value of 
0,—1.88 barns at v=2.2 km/s was used. The 
result was then 


2a, + a- 2 2.47 + 0.03) X 107? cm 
da 2L (0.73. 0.10) X 10-48 cm 


Only experimental points below 0.14 eV was taken 
into account. To determine the goodness of fit, 
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we also calculated pe For the two cases we got 
respectively 


X? = 29.9 Degrees of freedom = 25 


X = 33.1 Degrees of freedom = 26 


which is not too bad. 

In fig. 10 we show the total cross-section 
calculated from the last set of values and compared 
with the experimental data by Melkonian 
(1949). 


4.4 Oxygen. $ 
99,6 % of ordinary oxygen consists of the iso- 
tope O16, We have neglected effects of other iso- | 
topes. The O16-nucleus has spin zero and the elec- 2 
tronic wave function changes sign by inversion, 1 
therefore only odd values will be allowed for the 

rotational quantum number j. 

Curves for o,,, are shown in fig. 11 for j=11, 
and curves for o; for j—1, 11 and 21 in fig. 12. 
When averaging o; over initial states, to obtain the 
total scattering cross-section, j-values up to 25 
were taken into account. For higher neutron 
velocities, this average is to a good approximation 
equal to o; for j=11. 
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Fig. 8. The scattering cross-section for three initi, 
states for nitrogen gas, coherent transitions only include 
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The absorption cross-section is negligible for 


xygen. The magnetic moment of the O,-mole- 
ule gives, however, rise to a large magnetic scat- 
ering, which has been evaluated by Kleiner 
(1955). In fig. 13 a curve for the nuclear part 
of the scattering as calculated here is shown, and 
tlso a curve where the calculated magnetic scatter- 
ng is added and compared with experimental 
results by Palevsky and Eisenberg (1955) 
ind Melkonian (1949). 

To take benefit of the rather accurate data pre- 
sent in this case, we have subtracted the calculated 
scattering from the experimental values and 
compared the difference to the magnetic cross- 
section as calculated by Kleiner. As one sees 
from fig. 14, the agreement is very good. 

Our calculation of the nuclear cross-section also 
agrees well with the theoretical values presented 
for O9 by Halpern and Appleton (1953) 
for neutron wave-lengths above 5A. * _ 
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Fig. 9. Energy dependence of the coherent and the in- 
coherent scattering cross-section per N-atom for nitro- 
4 gen gas. 
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Fig. 10. The total cross-section for nitrogen gas. 
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Fig. 11. The cross-section for rotational transitions 


from the initial state j= 11 for oxygen gas. 
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Fig. 12. The scattering cross-section for three initial 
states for oxygen gas. 
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Fig. 13. The total cross-section for oxygen. 
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Fig. 14. The magnetic scattering cross-section as cal- 
culated by Kleiner compared with experimental data 
from which the nuclear scattering has been subtracted. 


5. Discussion. 


The validity of the Fermi pseudo-potential in 
conjunction with the Born approximation has been 
discussed by Zilsel et al. (1947) and by 
Lippman (1950). For hydrogen the cortec- 
tion will be of the order of 0.3 per cent. 

As can be seen from the graphs where we have 
compared the calculated total cross section with 
the experimental data (fig. 3, 10 and 13), the 
agreement is excellent for the lower and medium 
part of the energy region covered, but there is a 
deviation at the upper end. In the cases of Ny and 
Og this is certainly due to the neglection of effects 
from the zero point vibration. 

The deviation of 1 or 2 per cent for energies 
above 0.1 eV in the case of H, is, however, a bit 
difficicult to explain if the experimental results 
are assumed to be without systematic errors, since 
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e have corrected for vibrational effects using the 
irmonic oscillator model described above. 


Effects due to anharmonicity in the vibration 
ive been considered by Brimberg, but they 
ill be of second order and are certainly negli- 
ble when the parameters are chosen in the way 
e have done. 


For the high final rotational states, the centri- 
igal effect will be much more important. This 
ill cause a stretching of the molecule, thus redu- 
ng the overlap between final and initial wave 
inction and thereby the transition probability. On 
ıe other hand, the level spacing will be less and an 
creased number of final states is available for a 
rtain neutron energy. 


The effect of reduction will, however, be the 
ominating one, and a rough estimate based on 
irameters given by Herzberg (4950) shows 
at it would be of the order of 1 % at the upper 
id of the energy region considered hefe. „4 
Above 1 eV, one should also take into account 
gher order terms in the expansion of the spheri- 
1 Bessel function in eq. (49). 
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PHYSICA NORVEGICA 


Diffusion processes in the solid state at elevated 
temperatures have gained appreciable interest, and 
several methods for the study of such processes 
have been developed. Methods based on tracer 
technique are widely used, but various X-ray met- 
hods have also shown valuable results. Thus, for 
example, an X-ray absorption method was used 
by P. Lublin (1959) for the study of diffusion 
between aluminium and nickel. 

A special X-ray method, proposed by P. S. 
Rudman (1960), is based on the intensity dis- 
tribution of pairs of diffraction lines from a pow- 
dered, mixed, cold-pressed and heat-treated speci- 
men of a metal pair. A necessary condition for 
using this method is that the metals in question 
have the same crystal lattice structure and not too 
different lattice dimensions, so that complete reci- 
procal solubility is warranted. Each metal will, 
before heat-treatment, give sharp lines in a 
powder diffraction pattern. During the heat-treat- 
| ment, the metals will diffuse into one another and 
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Abstract 


A modification of a method, previously proposed by P. Rudman, for the study of inter- 
diffusion in metallic powder diffusion couples is described. The method is based on electrolytic 
deposition of one metal upon a foil of the other metal, and intensity distribution of X-ray diffr- 
action lines from the metal pair is measured by using the Guinier focussing principle, with a bent 
quartz crystal for transmission as monochromator. The setting of monochromator, specimen foil and 
counter is so arranged that the splitting of the K a-doublet is compensated. 


modify the lattice parameters, and this will in 
turn have an effect on the powder pattern, so that 
the pairs of lines, one from either metal, will gra- 
dually merge into one single line on prolonged 
heating at high temperature. The method has 
been applied for the study of interdiffusion be- 
tween copper and nickel by P. S. Rudman and 
B.Fisher (1961). 

There are apparently some practical difficulties 
and limitations connected with the method in 
this form. Thus, deformations and crystal strains 
may be introduced by cold-pressing in a varying 
and unknown measure, and the contact area 1s not 
known with any accuracy. Perhaps the most serious 
difficulty of the method, as described by Ru d- 
m a n, is that due to the influence of the splitting 
of the Ka-doublet on the line-width. Rudman 
has tried to eliminate this effect in two ways. One 
way is based on the application of a Fourier met- 
hod introduced by Stokes (1948), the other is 
a mathematical method proposed by Keating 
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Fig. 1. Guinier-de Wolff photographs showing the 111- 

and 200-lines for Cu-Ni interdiffused at various tempe- 

ratures for 4 hours. The diffraction patterns for unheated 

specimens are shown (above) for comparison in each 
case. 


(1959). Both these methods, however, are limited 
to the use of high order reflections. 

We have tried to modify R u d man’s method 
in several ways. First, we prepare our diffusion 
couples by starting from a foil of approximately 
0.01 mm thickness of one metal, on to which the 
other metal is deposited electrolytically. In this 
way an intimate contact between the metals is ob- 
tained, without introducing deformations or strains. 
Furthermore, the geometry of the diffusion couple 
becomes simple and unambiguous, and this is of 
importance for the subsequent mathematical treat- 
ment. Secondly, in obtaining the powder diffrac- 
tion pattern, we make use of the Guinier focus- 
sing principle with a bent crystal monochromator 
instead of the ordinary Debye-Scherrer method. We 
have, for a preliminary test, taken photographs of 
unheated and heat-treated specimens by the aid of 
an ordinary Guinier-de Wolff focussing camera. 
Photographs of some copper-nickel diffusion coup- 
les before and after various heat-treatments are 
reproduced in Fig. 1, from which it may be seen 
how the interdiffusion of the metals effects the 
intensity distribution of the diffraction lines, in 
this case the 111-and 200-lines. 

Very nearly the same principle is used for the 
accurate registration of the intensity distribution, 
but on a considerably larger scale. The principle 
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and arrangement is illustrated in Fig. 2. Rays from 
an X-ray tube are reflected from the radially 


oriented (1011)-planes of a curved quartz crystal 
plate, differing in this respect from the more 
usual Guinier-de Wolff camera. The radius of 
curvature of the crystal and the diameter of the 
Rowland circle is 400 mm. The crystal and the 
Rowland circle are turned through an angle equal 
to the Bragg-angle ©,, for the monochromator 
and radiation in question, so that the converging 
rays from the monochromator pass symmetrically 
about the center of the spectrograph circle, where 
the specimen foil is placed, as may be seen from 
Fig. 2. The monochromator crystal has a high resol- 
ving power, so that the Cu Ko-doublet is clearly 
resolved even though the Bragg-angle for Cu Ka- 
radiation is only 13,39. The half-width of the lines | 
is 2—3 minutes of arc, and the intensity amounts | 
to 2—3.10° cps. Diffracted rays from the specimen | 
are registered by means of a proportional counter 
fixed on an arm which is pivoted about the center 
of the spectrograph. The angular position of this 
arm may be adjusted by means of a micrometer 
screw and read with an accuracy of = 10 seconds 
of arc. The Bragg-angle can also be determined 
with a corresponding accuracy. 


Fig. 2. Arrangement for registration of the intensity dis- 
tribution for X-ray diffraction lines from metallic foil 
diffusion couples. The Bragg-angles Ba, Ayn. 0» and @ 
may be selected in such a way that the Ka-splitting is 
compensated (see text). يو‎ 


| The intensity is measured by the aid of a propor- 
[onal counter and a Siemens-Halske electronic 
»anel, equipped with a digital print-out. A counter 
lit of 0.3 mm is used in front of the counter in 
irder to obtain the necessary resolution. The num- 
ver of counts per minute is then reduced to 1000— 
£000, and it is necessary to count for several 
ninutes to obtain an acceptable statistical accu- 
acy. The angular resolution is very suitable for 
he low order reflection 111, the Bragg-angle for 
opper and nickel differing by 19 12’ (i.e. in 2 ©), 
hich corresponds to approximately 4 mm along 
he circumference of the spectrograph circle. A 
aigher angular resolution would reduce the inten- 
iities unduly, the intensity being spread over a 
arger area. The use of a monochromator has the 
idvantage of reducing the background conside- 
rably. a 

The method reported here offers a possibility 
pf preliminary investigations of the spécimen be- 
"ore heat-treatment, such as to crystal size distri- 
pution, dislocation density etc., by observing abso- 
ute intensity of reflection and the half-widths of 
the lines. The arrangement has furthermore one 
Histinct advantage as it removes the difficulty due 
٥٢ the splitting of the Ka-doublet in using high 
border reflections, which is inherent in the R u d- 
an technique. By using a monochromator cry- 
stal the rays are diffracted twice, and it is then pos- 
sible to obtain a compensation of the splitting. 
It may easily be seen from Fig. 2 that the con- 
Hition for such a compensation is approximately 
wiven by the relation: 


39, — Bragg-angle of Ka, — component from the 
metal pair, 


i — Bragg-angle of Ka, — component from the 
metal pair, 


— Bragg-angle of Ka, — component from the 
monochromator, 


mg 


— Bragg-angle of Ka, — component from the 
monochromator. , 


m1 


By proper choice of monochromator it is now 
bossible to select the reflection, which gives the 
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Fig. 3. Intensity distribution of 111-reflection with Cu 
Ka-radiation from (a) unheated Cu-Ni sample (b) Cu- 
Ni interdiffused at 750 °C for 1 hour. 


With Cu Ka-radiation one finds that equation 


(1) is very nearly satisfied for the reflection 1011 
from the quartz monochromator and the 111- 
reflection from both copper and nickel. This con- 
dition is verified experimentally by sharp and 
symmetrical intensity peaks for an unheated speci- 
men (see Fig. 3). Thus we are no longer limited 
to high order reflections, and we avoid the labo- 
rious mathematical methods of Stokes and 
Keating for correcting line-width because of 
the Ka-doublet splitting. 

The principle of this X-ray method is based on 
the fact that the two metals forming the diffusion 
couple, on interdiffusion develop a continuous 
distribution of lattice parameters. Consequently, 
each Bragg-reflection will have a form, of in- 
tensity versus angle, which depends on the con- 
centration distribution. Fig. 3 shows how the 
pure line profiles of the diffusion couple Cu-Ni 
(Fig. 3a) are changed on a heat-treatment for 
60 minutes at a temperature of 750 °C (Fig. 3b). 
Taking into account only those factors in the X-ray 
intensity expression that can vary significantly 
over the breadth of the pure diffraction line, one 
can find the X-ray equivalence of the so-called 


144 


concentration-penetration curves. From these cur- 
ves it is again possible to calculate how the diffu- 
sion coefficient varies with concentration and tem- 
perature. The calculation is based on Boltzmann's 
method for the solving the diffusion equation 
(Fick's 2. law), which leads to the following 
expression for the diffusion coefficient, D: 


Dene, =(1/2/)(dx/de),,[ * de (2) 


T 
where: 
; — time of anneal 
x — effective penetration 
c — concentration 


By means of equation (2) and the concentration- 
penetration curves (c-x-curves) it is now possible 
to find the concentration dependency of D by 
graphical integration. Equation (2) is strictly true 
only for infinite or semi-infinite media, and this 
means that the diffusion process must not be car- 
ried too far, in order that the metal pair can be 
regarded as infinitely thick. 

The results obtained by such a procedure repre- 
sent a bulk effect integrated over all directions in 
the crystals and their boundaries in the polycry- 
stalline specimen. The foil technique, described 
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here, is also suitable for investigating the detailed 
mechanism of the diffusion process by the aid of 
X-ray microradiographic methods. : 

The method has been used, for development of 
the technique, in the study of interdiffusion be- | 
tween copper and nickel in the temperature range 
from 600 °C—1050 °C, and for varying heating 
times. Similar studies on other systems, such as 
cobalt-copper and cobalt-nickel, are in the course | 
of development. A detailed report on this work 
will be published elsewhere. 
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Abstract 


An automatic apparatus based on a 4 x 4 coincidence unit for measuring the angular correla- 
tion of the annihilation radiation is described. The positronium formation is found to depend 
strongly on the crystallinity of the material and on the gas adsorbed on the surface. Thus NO, 
reduces the narrow component, H,, N, and argon show no difference compared to vacuum 
whereas O, by a conversion process increases the narrow component by 150 % at atmospheric 


pressure. 


Singlet positronium initially formed in the singlet state is slowed down to 0.5 eV and that 


formed by conversion to 0.1 eV. 


A method of data treatment by use of an electronic computer is described. 


. Introduction. 


The annihilation of an isolated positron-electron 
air is well understood on the basis of quantum 
slectrodynamics. 

The lifetimes of positrons in solids are long 
+nough to allow on the average almost complete 
lowing-down of a fast positron before annihila- 
ion. The electronic system of the solid will then 
fluence both the annihilation probability of the 
ositron and the spatial distribution of the anni- 
ilation quanta. Measurements of these quantities 
-an, therefore, yield information about the proper- 
lies of the solid. 

١ The two most important experimental methods 
ave been the study of positron lifetime in solids 
ind the study of angular correlation of 2-quantum 
nnihilation. 


*) Work sponsored jointly by Institutt for Atom- 
nergi, and Reactor Centrum Nederland, the Netherlands. 
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General reviews have been given by Simons 
(1958), Ferrell (1956) and Berko and 
Here old (1956). 


So far the study of positron annihilation has 
proved to be a tool of limited use in solid state 
physics, mainly because the interaction between 
the positron and the electrons of the solid is not 
fully understood. 

The purpose of this paper is to study the for- 
mation and dynamics of positronium (Ps, a bound 
state of a positron-electron pair in solids by the 
angular correlation method. 

As target material aluminium oxide was chosen 
because this compound can easily be obtained in 
different modifications. In the amorphous state it 
can also adsorb a large amount of gas, making it 
possible to demonstrate many features characteri- 
stic of a gas by means of a compact target. 

The principle of the experimental method is 
described in Sec. 2, and the apparatus used in 
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Sec. 3. The experimental results and a method of 
treating the data are presented in Sec. 4 and 
Sec. 5, respectively. A discussion of the results is 
found in Sec. 6, 


2. The angular correlation method. 
The general layout of the experimental setup is 
shown in fig. 1. Positrons from a radioactive 


SOURCE 
DETECTOR UNE دا‎ DETECTOR 
ee EE COE e 
i - NARGET B 
> < 


Fig. 1. Principle of experimental setup for measuring the 
angular cotrelation of annihilation quanta. 


source (64Cu or 22Na) incident on the target are 
slowed down and further annihilated. The target 
and the slits in front of the detectors will define 
an angle 0, and coincidences between the two 
detectors are recorded as a function of 0. 

If the linear momentum of the centre of mass 
of the electron-positron pair is p and that of the 
annihilation gammas k, and &,, respectively, we 
have for the case of a slowed-down positron (see 


Fig. 1). 


AL CE Yk, RI mc (1) 
The angle 0 is then given approximately by 
0 = p.lmc (2) 


where 5, is the component of p perpendicular to 
the plane through the two mutually parallel slits. 
There are two main cases: 


(1) The positron remains free. 


It is usually asumed that the positron is slowed 
down to thermal energies, i.e. a few hundredth of 
an eV before annihilation takes place. The anni- 
hilating electron will usually belong to the valence 
band and have a kinetic energy of a few eV. Hence 


P= Pp. +P- P- 
i.e. the angular distribution will be determined by 


the electron momenta. (For details see Stewart 
)1957(.( 
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(2) Positronium is formed. 


In 1/4 of the cases we get singlet (total spin 
zero) positronium. The most probable annihila- 
tion process is then selfannihilation (vacuum life- 
time 7,—1.25 x 10-10 sec) giving two photons. 


For the remaining 3/4 of the cases, in which 
triplet (spin 1) positronium is formed, direct 2- 
quantum annihilation is forbidden and we then 
have 3 competing processes. 


(a) 3-quantum annihilation, the vacuum life- 
time 7,—1.3 ۰ 10— sec, therefore this is 
an unlikely process in a solid. More prob- 
able is 


(b) «pick off» annihilation, ie. 2-quantum 
annihilation between the bound positron 
and a valence electron of the solid. If the ] 
solid contains unpaired electrons, Ferrell | 
(1958) has pointed out that 


(c) triplet Ps-atoms can be converted to singlet 
making the rapid self-annihilation possible. - 


The positronium atom is slowed down before 
annihilation takes place and will on an average 
have a momentum less than that of the valence 
electron. 

Self-annihilation of positronium therefore gives 
rise to a narrow component in the angular distri- 
bution, compared to the broad one originating 
from free positron annihilation and pick-off anni- 
hilation. 

The width of the angular distribution will be 
of the order of a few milliradians. Therefore a 
very good angular resolution is needed. 


3. Apparatus. 


A block diagram of the actual experimental 
setup is shown in fig. 2. Its basic parts have been 


=... 30 REGISTERS. 
16x30 SWITCH. 


+ 16 OUTPUT STAGES 
sree SCALE OF TWO ` 


~~ 


bnstructed and." described by G. Trumpy 
11960). There are 4 detectors on each side con- 
pu to a 4 x 4 coincidence circuit feeding 16 


utput channels. The arrangement is such that 
hese channels cover 16 successive values of 0, 
..6 mrad apart. 


The electronic circuits are standard and rather 
imple. To avoid a tedious calibration of their rela- 
've sensitivities, the coincidence rate for each of 
€ 30 angular positions which are observed is 
easured successively by all the channels, and 
ly the sum is recorded. This is achieved by mo- 
ing the counters on one side in steps equal to 
€ channel separation, 0,6 mrad. Simultaneously 
e switch between the output stages and the 
inal registers is shifted one step. The counting 
keriod in each position is 40 sec, apd one has to 
0 45 (= 30 + (16-1)) steps fortheand 45 
teps back, the cycling period will tMÉrefore be 
15 min, switching time included. 

This period is short enough to make any 
rther corrections for the decay of the source 
'64Cu, halflife 12.7 h) unnecessary, and errors 
jue to drift of the sensitivity are also greatly redu- 
d. 

The number of cycling periods required in one 
n can be preset and, once started, the operation 
s then fully automatic. 


The scale-of-twos are inserted (see fig. 2) be- 
use of the long deadtime of the mechanical 
gisters (0.1 sec). They are not read off. It can 
e shown that if N is the number of counts accu- 
mulated in a final register during a run and ø 
the corresponding standard deviation, then 
> o2/N > 1, the upper limit being approached 
vhen the chance for a count during the 40 sec 
eriod is much less than one and the lower limit in 
e case of many counts per period. 


The source is a copper foil of dimensions 
10 x 45 x 0.02 mm mounted on aluminium bac- 
ting. It is activated in the reactor to a strength of 
pprox. 0.1 curie and replaced every 24 hour. 

| The target is 0.1 g/cm? thick and 50 x 100 mm 
ide and is resting on a mica sheet. Both source 
nd target are placed horizontally with the long 
dge normal to the direction of observation. 

As indicated in fig. 2 the source cannot be seen 
rom the detectors, and the shield is also con- 
ed to prevent scattered radiation from rea- 


ing the detectors. 
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The lower part of the shielding is contained in 
a steel box (not shown in fig. 2) having a vacuum- 
tight joint to the upper part. The target area cap 
therefore be evacuated or filled with any gas of 
interest. 


4. Experimental results. 


In fig. 3 we show the angular distribution for 


COINCIDENCE RATE 


120.40 E87F-6E4 9724/02 226 85 ١ 
ANGLE /mrad 


Fig. 3. Experimental distributions for aluminium oxide 
targets subjected to heat treatment as indicated. 


three different targets of aluminium oxide sub- 
jected to various degrees of heat treatment as indi- 
cated. The raw materials was bóhmite, Al(OH)3. 
By heating under vacuum the water was rele- 
ased. X-ray examination showed that the target 
material corresponding to the upper curve was 
completely amorphous, for the lower one it was 
of rather good crystallinity, for the middle one it 
was somewhere in between. 
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VACUUM 


COINCIDENCE RATE 


TEA cox gam EE GEL ONE 
ANGLE / mrad 


Fig. 4. Experimental distributions for amorphous alu- 

minium oxide target with gas atmosphere as indicated. 

Dashed curve indicates broad component, fitted as des- 

cribed in text. Instrumental resolution curve shown on 
bottom. 


In figs. Á and 5 we show some angular distri- 
bution curves where we have varied the gas at- 
mosphere but kept the target material constant. 
These curves can be divided into 3 classes. 

_A. NOs. No narrow component seems to be pre- 
sent. The shape of the spectrum did not change 

. for pressure variations in the region 500—700 
mm Hg which indicates a "saturated" state. - 


ere is ‘evidently a narrow component present, 
and within the experimental accuracy it is the 
|» same in all four cases. AA MEM 


B. Vacuum, Hy, Ny and argon. Compared to A, 
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Fig. 5. Experimental distributions as a function ati 
O,-pressure. Broken curves fitted as described in text. 


5. Data treatment. 


The distributions presented in n B abov 
can be decomposed into three components, QU i 
background term, (2) a broad component. If th 
Ps-atoms were completely thermalized before anni 
hilation, this third term would be of a Ga 
shape. Also, if the thermalization is not com 
the shape should be close to a Gaussian if t 
is only one process active in the forn 
singlet Ps-atoms, and this we wi 
true for the vacuum او په‎ bu 


\ 
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| 
rq. (3) to the experimental points by the least 
iquare method. 

The function F(0) is fed into the computer in 
labular form, for intermediate values of 6, the 
ralue of the function is found by interpolation. 


For each experimental distribution, the para- 
meters A,, A», As l' and the zero point of 0 is 
calculated together with the corresponding stan- 
Hard deviations N C. given by eq. (4) and x? is 
also calculated. The distribution can also be prin- 
«ed out normalized with respect to Ay. 


The width of the narrow component, I, ob- 
rained by the procedure outlined above, will in- 
clude a contribution due to experimental broaden- 
ing. The angular resolution curve for the appa- 
ratus, as calculated from the geometry, is indicated 
in fig. 4. It is narrow compared to the observed 
curve and can be approximated by 3*Gaussian 
ith a width of ['„=0.53 mrad. As is welFknown, 
he folding of two Gaussian curves gives a new 
aussian curve, and the true width will be given 


y 


DEED: Js (5) 


. Discussion 


6.1 NO». 

As mentioned above, we assume the NOs- 
distribution shown in fig. 3 to contain only the 
broad component. Positronium is probably for- 
med in this case too, but is broken up again before 
lit is slowed down. Since the NO, molecule is 
known to be a very strong oxidizing agent, a sim- 
jple explanation is that the Ps-atom is oxidized, 
and a free positron is left. Another posibility 
is that either the positron or the whole Ps-atom 
forms a binding with the NO,-molecule. None 
‘of these cases would give a narrow compo- 
‘nent, but the shape of the broad one might be 
changed compared to the vacuum case. This has 
not been observed, however. 
| The broad component gives a good fit to the 
Gaussian curve exp(-0?|2 T3) with Dy = 4.60 
mrad. 


62 Vacuum, No, Ho and argon. 

Within the experimental error there is no dif- 
ference between the distributions in this case. 
When fitted to the function given by eq. (3) we 
| got 
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Ty, = (1:24 + 0 02) mrad (6) 
for the width of the narrow component. Correc- 
ting for experimental broadening by eq. (5), we 


get 


Leawiżmtad (7) 
If the momentum distribution is Maxwellian, the 
probability of finding a positronium atom with a 
momentum z-component between Pz and Pa t 


dp, is 


f») dp, = (4r mkT) * >< 


>< exp fi /4: 5 )8( 
where T is the corresponding equilibrium tem- 
perature. Since p,— 0 mc comparison with eq. (3) 
gives 

P= V 2kTime (9) 
Using the corrected value for T this gives for the 
mean kinetic energy 


3 
[c= z 27 = 0.48 V (10) 
This corresponds to a temperature of 3700° K and 
shows that the positronium atoms are certainly 
not thermalized in this case. 


The average value obtained for the narrow com- 
ponent fraction as defined by eq. (4) is 


N.C = 12 4 per cent (1.15 


Assuming that no conversion takes place and 
that the pick-off annihilation rate is negligible 
compared to the singlet self-annihilation rate, this 
means that close to 50 96 of the positrons form 
positronium. 

The fit of the theoretical function to the experi- 
mental data was good, but of course the number 
of points is too small to enable one to say any- 
thing definite about the shape of the narrow com- 
ponent except it is "Gaussian-like". 


6.3 Os and air. 

The enhancement of the narrow component in 
this case is no doubt due to triplet-singlet con- 
version by the two unpaired electrons of the 
O,-molecule. The effect has also been observed 
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by angular correlation measurements on gaseous 
O,-targets by He ta berg and Page (1957). 

To separate out the component due to conver- 
sion, we have again fitted the data to the function 
given by eq. (1), but now choosing the smoothed 
out vacuum distribution as the function F(0). As 
the average half width of the converted compo- 
nent we then got 


D*3(8 017557 002) tread (12) 


Correcting for experimental broadening as above, 
this becomes 


[° — 0.53 mrad (LS) 


The corresponding values of the kinetic energy 
and the temperature are 


E=011le/ T=830 K (14) 


Thus the energy of the converted positronium 
atoms seems to be a factor 4 less than of those 
formed directly in the singlet state. This result is 
of course reasonable, because the converted atoms 
have a total lifetime 7, 4- r, and the slowing down 
process can be more complete. 

In fig. 6 we have plotted the narrow component 


X X "TOTAL N. IC. 
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PER CENT 


200 400 600 800 
OXYGEN PRESSURE IN mm Hg - 


Fig. 6. The narrow component as a function 
of O,-pressure. 


fraction as a function of the oxygen pressure. The 
values for air at atmospheric pressure have also 
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been included, assuming that only the partial pres- 
sure of oxygen is important. 


The spread in the points is larger than the statis- 
tical errors, and is mostly due to difficulties in — 
getting reproducible results. The general trend of | 
the curves of fig. 6 is, however, thought to be : 
correct, and since quantitative results ate not of - 
much interest at this stage anyhow, the measure- 
ments have not been repeated. 
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It could be of some interest to extend the mea- 
surements to higher pressure, to see if the وا‎ 
molecule also exhibits oxidizing properties on the 
Ps-atom, but this cannot be done with our appa- 
ratus in its present form. 


What we have neglected in our way of trea- 
ting the oxygen data is that singlet positronium 
also is converted to triplet, thus the narrow com- 
ponent due to self-annihilation of positronium 
initially formed in the singlet state will be reduced. 
The results we have obtained for the width of the 
converted component should then be less than the 
true one. 


As a check, the Os-data were also fitted using 
the same function for the broad component as 
used in the vacuum case. The average value for the 
width of the total narow component at the hig- 
hest Og-pressure obtained in this way is 


43s 200 94 + 0.02) mrad (15) 


but the first method of treatment gave always a 
better fit. 

Comparison between eq. (6) and eq. (15) 
shows that the narrowing of the narrow component 
when O,-atoms are present is beyond doubt. 


7. Conclusion. 


We have demonstrated that the amount of posi- 
tronium formed depends greatly on the crystalli- 
nity of a solid target. This is also known from 
similar measurements on fused and crystalline 
quartz etc. : 


Our investigations on the role played by the gas 
atmosphere above the target indicate that, at least 
for the target material we have used, the positro- 
nium is formed in places where also the gas mole- 
cules can be present, i.e. positronium formation is 
a typical surface effect. : 


Lo 
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| ie other cases’ positronium may be formed in 
ores in the material, these need not be open to 
as molecules present at the surface. 


The explanation why the broad component is 
tot influenced by the gas atmosphere is probably 
hat the free positron is attracted by the negative 

) -ion of the oxide and not so much affected 
y the gas. For pick-off annihilation a polarization 
ffect of the same origin may be important. 

"The methods described here are therefore not 
nmediately applicable to for instance aquous solu- 
ions where one also introduces new ions. We 
Lm return to this question in a later paper. 
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Abstract 


The structure of the decay of the 14 keV level in 57Fe has been studied by delayed coincidence 
measurements for excited 57Fe nuclei bound in crystals of different thicknesses and temperatures. 
The data agree with earlier measurements, except for short times. 


Introduction. 3E M 27? (TM 
: R T 

When a nucleus bound in a crystal emits a y-ray, f exp (- 2L 0 E zT 5 | | 
there is in some cases a probability f different ry 
rom zero that the emission takes place without 
excitation of phonons in the lattice. This pheno- 
menon is known as the Mossbauer effect and by 
sing a Bravais model for the lattice and the 
IDebye theory one can show that 


Holland and his coworkers (Holland et al., 
1960 and Lynch eż al., 1960) have demonstra- 
ted the Mossbauer effect by observing the decay 
of the 14 keV level in 57Fe with separated absor- 
ber and emitter. We have performed a similar 
experiment with the absorber and emitter in one 


AD 1 k 3ER ( o ) | sample. 
parapi ET 
r Theory. 
ie > 2 It is a well known fact that the Mossbauer line 
PK szoł: i s from a thick emitter is considerably broader than 
ia د باو‎ oth © de, En — Faye the natural line-width T. Together with the un- 


certainty relation of energy and time, 7 PZ 


E, = energy of y-ray, M = mass of nucleus, 0 is Fourier transform of the Lorentz distribution 
the Debye temperature for the crystal, T the tem- prog 
N(EJdE= | (E — Ej)? چو‎ | dE gives the 
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decay spectrum N(1)dt—e r dt. In penetrating 

the sample the internal absorption will change the 

energy spectrum of the y-rays. As the ab- 

sorption coefficient of the y-rays is of the form 
2 


il 
A(E) œ [e-+] , we get most ab- 


sorption in the middle of the spectrum and the 
Lorentz curve will, by thick enough samples, be 
modified into a symmetric curve with a dip in the 
middle. Clearly this effect will change the time 
distribution, and is has been calculated in an app- 
roximate quantum mechanical manner by Harris 
(1961) and by means of a classical model by 
Hamermesh (1960) and Lynch ef al. 
(1960). (The most important approximation is to 
use a plane wave representation for the y-ray.) 

The result for one single spectral line with 
Mossbauer coefficient f=1 is 


la) oe] Ju (Bt X37) | um an 
where | a(t)|2 is the time-dependent transmission, 
A-1[r, and f is a thickness parameter given in 
Sec. 4. 
When we take into account the hyperfine struc- 
ture and the relative intensity of the 57Fe lines, 
and that 0 < f « 1, we get |a()|* — اا‎ 


x | 1 -fsh DU 2 2F اف‎ la 
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Fig. 1. Decay of 6 


3. Experimental Arrangement and Procedure. 


Measurements were made to show how the in- 
ternal absorption changes the decay of the 14 keV 
level in 57Fe by measuring the time correlation be- 
tween the 122 keV and 14 keV y-rays from the 
decay of an excited "'Fe nucleus (Fig. 1) in a 
"Fe lattice. The experimental arrangement is 
shown in Figs. 2 and 3. The samples were made 
from polished soft iron, on which a solution of 
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Fig. 3. Experimental arrangement. 


57Co Cl; was evaporated at room temperature 
and thereafter heated to 900*C for two hours in 
a thin hydrogen atmosphere. To avoid surface 
contamination the HCI attached surface was re- 
polished, and a layer of 57Fe enriched to 78 96 
was electrolyzed onto the sample. The layers were - 
2.2 and 5.0 mg 57Fe/cm2, and measurements 
were made at room temperature for both samples - 
and for the thickest sample at liquid nitrogen 
temperature also. The counters for the 122 and 
the 14 keV y-rays were Nal crystals, 0.08" and 
0,01" thick, mounted on an RCA 6342 and a 
56 AVP photomultiplier, respectively. The time- 
to-amplitude converter assembly was the same 
(except for small modifications) as previously 
used for time-of-flight measurements of fission 
neutrons (Skarsvág and Bergheim, 1963). 
The integral discriminator setting for the 122 ٧ 
y-tay was put well above the 14 keV peak, and 
the single channel analyzer was set to accept 
the photopeak of the 14 keV y-ray. The zero time 
was determined either with a 60Co source, or by 
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iscriminating the 14 keV y-ray and using a bad 
ometry to get radiation scattered from one 
punter to the other. From the full width at half 
aaximum of the peak, and from the steep rise 
t the left side of the exponential curve for pure 
Co, the time resolution was found to be about 
nsec. 

The calibration of the time scale was checked 
very day. An overall check of the equipment was 
aade by measuring the exponential decay of the 
keV level of 57Fe (Fig. 4). The line corres- 
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Fig. 4. Decay of the 14 keV level ۶ه‎ 6 


»onds to the half-life of 10-7 sec. The background 
vas determined by using different 57Co sources 
r each counter. The counting rates in the back- 
tround measurements were adjusted to about 

e same as in the actual measurements. The back- 
;round was about 2 % for the 2,2 mg/cm? and 
10 96 for the 5.0 mg/cm? sample. 


. Results. 

١ The number of counts in each time channel, 
nultiplied by e, is plotred in Figs. 5—7, 
where the errors are standard deviations. The 
lashed curves are theoretical with a thickness 
»arameter B — No, f, which is the expected value, 
here N is the number of 57Fe nuclei per cm?, 
7, 1.48 - 10-18 cm2 is the peak absorption cross 
ection for no hyperfine splitting and f—0.6 at 
oom temperature and f=0.7 at liquid nitrogen 
mperature. The fulldrawn curves are for" oS 
No, f which is the value used by Lynch e/ al. 
1960). 


Internal Mossbauer-effect in 57Fe 


200 nsec 
TIME 
Fig. 5. Time spectrum (with ordinate multiplied by e?7) 
obtained with absorber of thickness 2.2 mg ®*Fe/cm? 
at room temperature. The solid curve is for 8 = 2N oo f 
and the dashed curve for 8 = N o, f. 


1 2 E 
00 TIME 00 nsec 


Fig. 6. Time spectrum obtained with absorber of thick- 
ness 5.0 mg 9'Fe/cm? at room temperature. The solid 
curve is for 8 = 2N cpf and the dashed curve for 6 = 
No, f. 


We think that the internal effect could have 
some interest in testing whether the Mossbauer 
effect arises when the lines are extremely narrow 
or whether line-broadening effects such as the 
spin-spin interaction make the effect disappear. 
In an ordinary Mossbauer experiment several 
effects can make the resonance vanish. By per- 
forming experiments in one sample, which we 
have shown is feasible, it should be possible to 
avoid the unwanted energy shifts to a very large 
extent. 
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100 
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Biz. 7. Time spectrum obtained with absorber of thick- 
ness 5.0 mg 57Fe/cm? at liquid nitrogen temperature. The 
solid curve is for 8 = 2N 6) f and the dashed curve for 
B-—No,f. 
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Abstract 


Absorption coefficients for pure copper are measured for wave lengths slightly above and below 
the CuK-absorption edge. Measurements of absorption coefficients, relative to those for the CuKf,- 
line, are made for the peak of the CuKf2,;-band in pure copper, some copper alloys and copper 
compounds. It is shown that the absorption of the CukKf » ; -line exceeds the normal value in pure 
copper and in copper alloys and may vary with the state of alloying. The absorption for the 
CuKf»,s -line is normal or only slightly above the normal value in the case of copper compounds 
as absorbators. The variation of the absorption of the CuKf»,;-line in copper alloys is connected 
with a displacement of the absorption edge and of the position of the Fermi level of copper. 


1. Introduction. 


The study of the X-ray emission bands and the 
corresponding absorption spectra from metals and 
alloys is a way of obtaining informations as to 
the distribution of energy states for the outer elec- 
trons in such metals and their alloys. The short 
wave length emission band of an X-ray series 
originates from electron transitions from 6 
occupied part of the outer energy band or over- 
lapping bands to an inner ionized level. The 
distribution of intensity versus frequency, of 
versus photonenergy, in the emission band there- 
fore represents the electron density distribution in 
the occupied part of the outer band or overlapping 
bands, only modified by the transition probabili- 
ties and the rather small width of the inner 


` ionized level. The absorption edge and its nearby 
` structure, on the other hand, correspond to elec- 


tron transitions from an inner level to the first 
unoccupied energy levels satisfying the selection 


rules. Thus, the absorption spectrum reflects the 
electron energy state density distribution of the 
unoccupied part of the outer energy band. 


In the case of the elements in the range of 
atomic numbers from about 20 to 35 the short 
wave length emission band in the K-series will 
be either KB, or Kf, corresponding to transitions 
K—Nyrni or K=Mqypy respectively, and the 
transitions on absorption can be either Nr; gi —K 
or Mi, y —K, depending on the atomic number 
and thus the number of electrons in the outer 
shells. 

The difference in wavelength between the 
absorption edge and the emission band peak is 
found to be very small as long as the outer shell 
in question is only partially occupied by electrons. 
This wave length difference must, however, be 
expected to increase rather abruptly with atomic 
number from the element where the outer shell 
for the emission transition for the first time is 

\ 
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fully occupied, as the electron on absorption in 
that.case must be raised to the next higher shell. 

In the case of copper it is not possible to dis- 
criminate between the 8, and the 9;-band because 
of the overlapping of the outer shells ۰۷ 
and ور‎ The short wave length emission band 
is therefore designated as the Kf, „-band. This 
band has a sharp intensity peak which is here for 
brevity called the K5, s-line. 

The CuK-emission and absorption spectra have 
been investigated by several authors, first by 
Lindh (1930), Veldkam p (1933) and later 
by Swedenborg and Claeson (1937) and 
Kurylenko (1940). Very accurate measure- 
ments on the absorption spectra of the elements 
from Z=26 to Z—35 were carried out by Bee- 
man and Friedman (1939). Their results 
are, however, given in the form of curves on an 
arbitrary scale for log I/I, versus wave length, 
and absorption coefficients for the CuK$, ;-line 
cannot be found accurately from their measure- 
ments. Very extensive works on the K-absorption 
edge and the K£,—K$8, „-emission spectra of 
copper were carried out more recently by Eda- 
moto (1950) and by Hayashi (1936, 1952). 
In all these works accurate data for wave lengths 
and details of the spectra are found, but accurate 
values of the magnitude of the absorption coeffi- 
cients for the CuKf, ;-line are not given. 

The absorption edge is known to be displaced 
slightly for an element as constituent of an alloy 
or of a chemical compound compared with that 
of the pure element. Thus the displacement of tke 
K-edge is measureable for elements of low atomic 
number with Z «30. It is, however, thought 
that it would be easier to observe a change in the 
absorption coefficient of a fixed line very close 
to the edge, than to observe directly a displace- 
ment of a more or less complex absorption spect- 
rum because of alloying or/chemical binding. 
Provided that the form of the absorption edge 
can be approximated analytirally, it will then be 
possible to obtain the displacement of the edge 
from a measured change in the absorption for an 
emission line lying very close to the edge. 

In this paper a preliminary report is given on 
measurements of the absbrption in pure copper 
of wave lengths around its K-absorption edge 
and on measurements of the absorption of the 
CuK, s-line in pure copper and in some copper 
alloys and chemical compounds. 
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2. Experimental data. 


A sealed off Philips X-ray tube with a copper 
target was used as radiation source in this measure- 
ments. The take-off angle of the radiation from 
the target surface was approximately 6°. This 
take-off angle may be of some importance here, 
because the maximum of the Kf, ;-emission band 
might be slightly dependent on this angle due 
to self-absorption in the target material. The tube 
was run under ordinary condition of 35-40 kV _ 
and 12-15 mA. | 

A spectrograph, built in this Institute (So rum, _ 
1947), and equipped with a quartz concave crystal 
for transmission (Cauchois, 1935), has been 
used in this measurements. The diameter of the 
Rowland circle is 400 mm. The crystal plate is 
inclined to the Rowland circle in such a way so 
as to take advantage of the parallel focussing 
principle developed by Fankuchen (1937). 
Measurements were made using both 101, 102 
and 203-reflections from the quartz crystal, having 
relative intensities of 100 : 12:11. The 203-re- 
flection is supposed to give the most accurate 
results because of its higher resolution. Most of 
the results reported here are therefore obtained 
by using this reflection. The resolution of the 
monochromator is quite good, thus, f.ex. the 
CuKa; ,-doublet is clearly resolved even for the 
101-reflection, having a Bragg-angle of 13,39. 
The crystal plate is cut in such a way that the 
inclination between the tangent of the Rowland 
circle and the crystal plate is 47?, 26? and 199 
for the reflections 101, 102 and 203 respectively. — 

The X-ray intensities were recorded by the aid 
of a xenon-filled proportional counter and a 
Siemens-Halske electronic panel, using a pulse 
height analyzer with channel width of 1 V. Counts 
per minute were printed automatically by the aid 
of a Kienzle digital print-out, and amounts to 
10000—20 000 for the maximum of the Kia 
band without the absorption foil being inserted. - 
The slit in front of the counter was 0,09 mm - 
wide, corresponding to a width of 20 second of — 
arc in the angle 6. The 26-angle was adjusted by 
means of a micrometer screw arrangement and 
could be read from a graded scale with a vernier, 
which allowed readings to be made with an 
accuracy of —10"-— — LT 

The thicknesses of the absorption foils were, 
according to a suggestion by Nordfors (1960), 
chosen so as to give optimum conditions for the 
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E ت‎ 45 of' absorption coefficients as far as 
ossible. The thicknesses of the foils were deter- 
mined by weight, or by measurement of the ab- 
sorption of the Kei- and/or the Ka, ,-lines in 
rases where absolute determinations of the coeffi- 
cients were made. The foil was placed between 
cylindrical surfaces of the same curvature as the 
crystal holder in such a way that the X-rays every- 
where pass perpendicularly through the foil. 


The 8,-line and the B5;-band in the K-emission 
spectrum from the copper target are registered 
y the aid of the curved quartz crystal with the 
03-reflection in transmission. The width at half 
maximum intensity of the 8,-line is in angular 
measure 2,2. The angular difference between the 
eflections of the B,-line and the Bs line is, 20° 


„w „5 
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in 2 0, There is a slight broadening of the lines 
at the bottom, and this is presumeably mainly due 
to air scattering. 


3. The CuK-absorption edge and 
its nearby structure. 


The absorption edge of pure copper was 
registered with a copper foil of thickness 
35,94.10-4 cm., as determined by weight, and 
checked by measuring the absorption of the 
CuKa, „-lines. In some cases the measurements 
were also carried out with a foil of twice this 
thickness. The absorption coefficients (total mass 
absorption, u/p versus wave length A) were deter- 
mined by using the spectrum of the CuK8, and 
the CuK, as primary radiation. The obtained 
absorption coefficients are given in Table 1. 


d Table 1. 

bserved and calculated mass absorption coefficients of pure copper for wave lengths around the K—absorption edge 
of copper. v is frecuency, R the Rydberg constant, u/p mass absorption coefficient, T is the halfwidth of the K—level, 
ke is the effective halfwidth corrected for experimental broadening. 


Wave length 4 | 


A»/R 2n E (eV) 5 


Mass absorption coefficients 4/0 


: : v|R from abs. from abs. calculated | calculated 
enger. | edge | edge | observed with Tz/2 = | with 122/2 = 
| 0,725 eV 0,9 eV 

1,40367 649,20 —-11,06 —149,8 38,6 

1,329710 652,26 — 8,00 —108,4 38,2 

1,39492 653,28 — 6,98 — 94 6 38,2 

133 92278 ,) 654,57 — 5,69 — 77,1 S MAS) 

1,39042 655,39 — 487 — 660 38,1 37,84 38,16 
1,38824 656,42 — 3,84 — 52,0 38,6 38,11 38,40 
1,38618 657.40 — 2,86 — 38,74 38,8 38,54 39,00 
1,38503 657,94 — 2,32 — 31,43 39.9 38,90 39,48 
1,38400 658,43 — 1,83 = ANTO 40,7 39,45 40,18 
1,8290 658,96 — 1,30 == WIG 41,4 40,48 41,47 
115259222 659,28 — 0,98 — 13,28 42,2 41,63 42,28 
1,38181 659,47 — 0,79 — 10,71 45,8 43,07 44,43 
1538157. 659,68 — 0,58 —. 7,86 51,0 45,19 47,46 
1,38094 659,89 — 0,37 — 5,01 DOS. 50,07 32,95 
1,38050 660,10 E O16 — 0 TLE 65,15 71,06 
1,38006 660,31 + 5 + 0,68 196 239,2 230,0 
1,37961 660,52 -|--0;,26 + 3,52 243 292,9 288,4 
137919 660,73 E OA A „6:34 185 300,2 297,9 
1,37875 660,94 -- 0,68 + 9,21 303 

1,378231 661,15 + 0,89 206 167 

E. 1,37744 661,57 Ba n WITS 235 
i 1,37307 663,67 E 221 + 46,20 203 
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Several corrections have to be applied in cal- 
culating the absorption coefficients from the ex- 
perimental observations. The applied corrections 
are: 

i) Air-scattering of the f,-line. It is found 
that the air-scattering varies approximately with 
the third power of the angular displacement in 
minutes of arc from reflection position of the 
B,-peak, and it is proportional to the intensity of 
the B,-line. The corresponding correction in the 
intensity passing through an absorption foil is 
calculated by using the absorption coefficient of 
the B,-line. 

ii) A similar correction is applied for air- 
scattering of the 8, s-line. As the intensity of this 
line is very weak this correction is very small 
except for very close to the 8, ;-peak. 


iii) The background of continuous radiation 
is rather low as the pulse height analyzer will 
cut away second and higher orders of spectra. 
The background is found to decrease approxi- 
mately linearly with increasing angular position 
over the small range of angle which has to be 
measured. The magnitude of this correction is 
found from the measurements on the short wave 
length side of the absorption edge. 


iv) Correction due to background from cosmic 
rays, stray scattered rays etc. has been applied, 
but this correction is very small, usually below 
10 cpm. 

v) In the case of an absorption foil containing 
nickel (see following section) it is found necces- 
sary to apply a correction for secondary radiation 
from Ni, because of the strong absorption of the 
B-lines in Ni. An approximate value for this 
correction is found from the deviation from 
normal absorption slightly above and/or slightly 
below the f-line. 


No correction has been applied for scattering 
of the primary X-rays in the absorption foil. The 
results obtained here therefore represent the total 
absorption. 

The angular displacement (in 26) between the 
peak of the K£, „-band and the inflection point 
of the absorption edge is 1,72 minutes of arc, 
corresponding to an energy difference of 5,0 eV. 
This value is comparable with the value of 4 eV 
as given by Blochin (1957). The main width 
of the absorption edge is to some extent deter- 
mined by the experimental broadening, parti- 
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cularly by the slit width. There is, for the same 
reason, only a slight indication of the absorption 
line on the long wave length side of the edge, 
as observed by Beeman and Friedman Ics 
The absorption is found to vary strongly and 
rather irregularly on the short wave length side 
of the edge. This must be mainly due to the fine 
uructure of the edge, but as the absorption here 
is very strong and the number of counts low, 
statistical fluctuations may also be responsible for 
this variation to some extent. It is to be noted 
that the scope of this work is to measure the 
absorption of the Kf, s-line, and the foils used 
are therefore too thick to register the fine struc- 
ture accurately. 


4. Absorption. of the CuK/.; - line in pure copper 
and in some copper alloys. 

In order to observe a possible displacement of 
the absorption edge on alloying, the absorption 
for the peak of the Kf, „-band is measured for 
pure copper and for certain copper alloys. These 
measurements were carried out in such a way that 
the arm carrying the crystal monochromator and 
the arm carrying the slit and the proportional 
counter are left in fixed positions, while the 
absorption foils were inserted successively for 
measurements. This is to acertain that all measure- 
ments are taken at exactly the same angular position 
and thus for exactly the same wave length of the 
radiation. Measurements for pure nickel and for 
pure aluminium were also done merely as a control. 

The calculation of the mass absorption coeffi- 
cient p/p or, in cases were the thickness (x) is | 
not known, the value of (u.x), is carried out in 
the same way as described in the -preceeding 
section. 

The obtained results are set out in Table 2. In - 
the case of an alloy the total linear absorption — 
coefficient for the copper content has been cal- 
culated from the knowledge of chemical compo- - 
sition. The percentage increase in absorption of 
the KB, s-line above the «normal» absorption is | 
given for comparison in the last column. This 
anomality is further discussed below. In order to - 
bring out the effect of alloying more clearly a 
similar series of measurements was carried out 
for a wave length about halfway between the — 
KB, -band peak and the inflection point of the 
edge. The results are given in Table 3, and it 
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Table 2. 


kbsorption of the CuK/f2,s—line (A v/R& —0,4 on the short wave length side of the edge) in pure copper and in some 
pper alloys, measured relatively to the absorption of the CuKf,—line. 4 is the absorption coefficient and x is the 
hickness of the specimen. 


pa (ux)a Coefficient Percentage 
for the for the a fi i 
Absorbator j " Ha fOr increase 
CuK p\— CuK د‎ al ey CuK 5و‎ above «normal» 
line line in Cu absorption 


1,2138 1,6419 1,353 445,2 36,6 
1,2591 1,5585 1,238 473,8 42,9 
rass (85Cu15Zn) 0,8719 1,1635 1,335 475,0 43,3 
| - bronze (90Cu10AL 0,8820 1,2084 1,371 463,5 39.7 
ynite (95,5 Al 4,5 Cu) 0,8512 0,8696 1,022 466,0 40,5 
i— Cu foil (unheated) 3,4744 3,6800 1,059 
«  (4hat 600°C 3,6475 3,8072 1,044 
* fab at 67 3,5152 3,6960 1,051 
« (4hat 850°C 3,1003 3,1818 


« (3h at 1050°C 2,5982 «| = 2,8000 
Aluminium 1,2069% — 1,1857 
| Nickel 3,2327 3,2000 


may be seen that the absorption is considerably not be perfectly even, and this may influence the 

increased for the alloys as compared with that of results to a certain extent. The results, which are 

pure copper. set out in Table 4, clearly brings out, however, 

that the absorption of the K£, s-line in such 

compounds is nearly normal or only slightly above 

5. Absorption of the CuKéf2s—line in the normal value. The displacement of the X-ray 

copper compounds. absorption edge must therefore go in the opposite 

Similar measurements have been performed for direction as compared with that observed for the 

some compounds of copper. The absorption alloys. This result for chemical compounds are in 

«foils» were in this case prepared from powdered agreement with those earlier found by Sanner 

samples, which were pressed between thin collo- — (1941) and others from direct observation of the 
dion foils. The thicknesses of these specimens may displacement of the edge. 


. Table 3. 


Absorption of the CuKf»,;—line (Av/R& —0.25 from the edge of copper) in pure copper and in some copper alloys. 
u, and us are the linear absorption coefficients in cm-! of the CuK2, and the CuK 2,s—line respectively. X is the 


| thickness of the absorbator. u, and(u/o), are calculated for the copper content. 


Percentage 


increase 
Absorbator eee AA mt (ulọ) » for Cu | above «normal 
value 
| 1252 6) els واه د نه‎ aa ee E akt TERA pape, 
| Copper 1,2138 1,8412 511 57,5 56,8 
Brass (67Cu33Zn) 1,2591 1,7159 546 61,2 67,5 
` Brass (85Cul15Zn) 0,8719 1,3145 547 61,4 67,8 
` Al—bronze (90Cu10AI) 0,8820 1,3928 20 60,0 61,8 
Lynite (95,5 Al 4,5 Cu) 0,8512 0,8848 45224 58,5 60,2 
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6. Discussion. 

a. The absorption edge. The form of a true 
absorption edge was investigated theoretically by 
Richtmyer, Barnes and Ramberg 
(1934), who predicted an arctangent curve for 
the variation of the absorption coefficient across 
the edge and they deduced the expression: 


T TCO r | (E. ده‎ Eon] xis $ | | 


whete Too is the absorption coefficient for the 


(1) 
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widths of inner levels could be obtained from the 
half-width of absorption edges, provided that 
these are true absorption edges. Some measure- 
ments of the width of the K-level of various 
elements have recently been made by Meisel 
and Nefedow (1961), from which work a 
value of 1,45 eV is obtained for the K-level of 
copper. This value of T; has been used for the 
calculations given in the sixth column of Table 1. 


The relation deduced by Richtmyer et al. 
and given by Blochin is, however, to be 


Table 4. 


Absorption of the CuK»;—line (Av/R=—0,4 on the short 


compounds measured relatively to the absorption of the CuK, : | 
B2s—line respectively. X is the thickness of the specimen. (u) and (u/e) are calculated 


of the CuK$8, — and the CuK 
for the copper content. 


wave length side of the absorption edge) in some copper 
—line. u, and u, are the linear absorption coefficients 


Percentage 
(ux), for (ux) for Ua (cm-1) for | (u/g)s for the | increase above 

Absorbator ‘Cakes, CuK$»;5 the Cu content Cu content «normal» 

absorption 
CuO 1,944 15002 343,4 38,5 5,4 
O Cl 3,820 3,939 366,0 41,0 1273 
Cuci, 22771 3,048 346,4 38,8 6,3 
CaCO; 1,830 21071 382,3 42,7 14,8 
CuSO,.5H,O 0,424 0,449 395,0 44,3 18,7 
(CH OOO) CTH O 0,574 0,613 374,0 41,9 12:3 
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short wave length side and away from the edge, 
Eiofo is the photonenergy for transition of an 
electron from the distribution maximum of the 
inner level to the limit of the continuous energy 
region, or in other words, to the Fermi level. 
Ex is the variable photonenergy and finally T; 
is the width of the inner level in question. The 
direct application of this relation to absorption 
edges is, however, in many cases somewhat 
questionable because the edge may be, as shown 
by Beeman and Friedman 1.c., complicated 
by the superposition of more or less pronounced 
absorption lines. This is also the case for copper. 
The experimental broadening is so marked, how- 
ever, that the faint absorption line is hardly 
visible. The edge, therefore, approximates fairly 
close to an arctangent curve. The value of Ty 
for the K-level of copper is needed for the use 
of equation (1). Accurate determinations of the 
widths of inner levels are difficult to obtain and 
rather scarce in the literature. In principle such 


changed slightly for this calculations in order 
to give correct values on both sides of the edge, 
thus: 


tan"! 


1 


Esta] 


It is then easily seen that 7, = rę for E, >> 
وون‎ te~ ry for B€ < En and Tr.=$ 
(7 oo + 79) for E, - Ej, n. The values of the absorp- 
tion coefficient can be taken as 37,5 and 310 for 
the long wave length and the short wave length 
side of the edge respectively. r is here the total 
mass absorption coefficient (uo). 

It may be seen by comparison of the calculated - 
values of the absorption coefficients in column 6 
of Table 1 with the observed values in column 53 
that the agreement is fairly good 30—70 eV - 
below the edge but the calculated values deviates - 


= 


hore and more from the observed ones as the 
dge is approached. This might be due to a too 
pw value for the width of the K-level (T;), 
jut is more likely to be caused by the experimental 
itoadening. This experimental broadening may be 
aken into account by adding the experimental 
ialf-width, which is mainly determined by the 
lit width and lies between 1,0 and 1,2 eV, to the 
width of the inner level T; in the relation above. 
this approximation the effective width of the 
vel (Tye) is thus taken as T,2—T;-F T... 
ind equal to 1,8 eV. Absorption coefficients cal- 
lated with this value of Ty, are listed in the 
eventh column of Table 1. It is obvious that these 
ialues are in better agreement with the observed 
ines than those i the preceeding column. There 
still some discrepancy for one or two values 
llose to the edge, but this might possibly be 
aused by the mentioned absorption;*line for 
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b. Tbe influence of alloying and of chemical 
ending. Assuming that the third power law for 
e absorption coefficient versus wave length is 
ralid in this wave length range, the coefficient 
ior the CuK, ;-line would be 2,4 % lower than 
ior the K8,-line. It is then clear from the present 
»bservations that the absorption coefficient of the 
Bə s-line in copper and copper alloys exceeds the 
«normal» value to a considerable extent, thus for 
bure copper and copper alloys by 30—40 %. In 
'opper compounds, on the other hand, the absorp- 
ion is nearly normal or exceeds the normal value 
pnly by a few percent. It is also seen from the 
measurements with pure copper as absorbator that 
1e deviation from «normal» absorption extends to 
L wave length lying about half-way between that 
of the Bı- and the f, s-line, though the deviation 
lhus far from the edge becomes very small. 


= 


| The results obtained with Cu-Ni-foil as absor- 
mator clearly shows that the absorption of the 
JuK$8, s-line in copper depends on the state of 
lloying. The foils are here the same in every case, 
but the interdiffusion of copper and nickel, and 
thus the state of alloying, varies with temperature 
nd time of anneal. The Cu-Ni-foils are, however, 
father unfavourable for measuring the absorption 
wery accurately, as the CuKf-lines are very strongly 
libsorbed in nickel. The absorption coefficients for 
pure Ni and Al are, as expected, lower for the 
CuK, s-line than for the CuKf,-line. 
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In the case of copper in brass an increase in 
the absorption coefficient for the CuKf, ,-line 


is observed compared with that of pure copper. 
Such an increase is also observed for the Cu-Al 
alloys. 


Equation (2) may be solved for Es—Eiofo = AE, 
in order to calculate the displacement of the 
absorption edge from observed values of the coeffi- 
cient for the CuK8,,-line in the specimen in 
question. Thus we get: 


— T/2 — 7 ool 2) 
TT 


A Gr 
AE, = (Trel 2) tan | T (3) 


OOM 0 

When the value of AE, for pure copper is sub- 
tracted from the value for an alloy or a compound 
we get the displacement of the edge relative to 
that for pute copper. The results of the calculati- 
ons are set out in Table 5, where the displacements 
of the edges have been converted to wave length 
differences. The displacements of the edge for 
chemical compounds are in the direction of 
shorter wave lengths, the magnitude of the dis- 
placements observed here, are, however, somewhat 
larger than those observed by Sanner E. GO. 
especially for CuO. This may possible be due to 
an influence of the said absorption line for copper. 


The displacements of the edge for alloys are 
smaller than for compounds and lie in the direc- 
tion of longer wave lengths. The displacements 
for both specimens of brass are the same within 
the limits of accuracy, and amount to — 0,7 eV. 
For the Cu-Al alloys the displacements of the edge 
go in the same direction, but are only about the 
half in magnitude. Both Cu-Al alloys are in fact 
two-phase alloys and copper must be assumed to 
be saturated with Al in both cases. 


A displacement of the absorption edge of an 
element on alloying could possibly be correlated 
with a change of the Fermi level of the same 
element. Such a correlation is, however, not à 
simple one, especially since the Fermi surface of 
the noble metals are shown by Pippard (1957) 
and others to be far from spherical. According to 
Cohen and Heine (1958) there are two 
effects upon alloying. Firstly, there is the filling 
up of the band, and secondly, the distortion of the 
Fermi surface will be affected. The distortion may 
be increased or decreased. In copper it is most 
likely to be decreased upon alloying with zinc. 
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Table 5. 


Displacement of the absorption edge due to chemical binding or alloying of copper (in A units), and change of the 


Fermi level on alloying of copper in eV. 


Mass absorption : Change of Fermi Change of Fermi 
coefficient in Displacement level from meas. level from meas. 
Specimen the copper con- of edge m A at APIR —0,4 at Av/R24—0,25 
E dem wu from the edge from the edge 
CuO 38,5 — 0,0067 
CuCl. 41,1 —0,0021 
CuCl, 38,8 —0,0055 
Cus GO; 42,7 — 0,0013 
CuSO,5H,0 44,3 —0,0008 
(CH COO), Cu. HO 41,9 — 0,0016 
Brass (67Cu33Zn) 53,1 +0,00010 —0,71 
Brass (85Cu15Zn) 22:2 +0,00011 —0,68 
Al—bronze (90Cu10A1) 52,0 4- 0,00005 —0,41 
Lynite (95,5 Al 4,5 Cu) W - 0,00006 


The Fermi level of copper has been determined 
in several ways, thus, f.ex. a value of 7,4 eV 
was obtained from the width of the La-emission 
band (Sórum, 1947). On the assumption that 
the displacement of the absorption edge is mainly 
due to a change in the filling up of the band, 
the Fermi level would be reduced to 6,7 eV in 
a-brass, in spite of the apparent increase in the 
number of conduction electrons on alloying with 
zinc. This is somewhat surprising, although the 
suggested bound states for the solute atoms 
(Friedel, 1954) and the scheme for metallic 
valences put forward by Pauling (1948) 
would indicate that zinc atoms do not contribute 
by as much as two electrons per atom to the 
conduction band. It may be however, that a change 
in the distortion of the Fermi surface, a dilation 
of the lattice and possibly a change in the structure 
of the absorption edge are also involved. 
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with the theory of Bardeen e al. (1959). 


conductivity above T, was established. 


. Introduction. 


Both electrons and phonons contribute to the 
ermal conductivity in metals, and interpretation 
»f the observations is complicated because the two 
ontributions can be limited by several different 
icattering processes. 

In a superconductor near the transition tempe- 
ature (Te) the main part of the heat will be 
karried by the electrons. Below T, the electronic 
bart of the thermal conductivity will decrease 
xhile, on the other hand, the lattice contribution 
will increase. This can be understood qualitatively 
‘rom a two-fluid model assuming that only normal 
-Jectrons scatter phonons and that this number 
lecreases below Te. The initial rise in the lattice 
bart of the thermal conductivity will result from 
Ihe decreasing numer of electron scatterers for the 
phonons. 
| Based on the Bardeen, Cooper and 
Schrieffer (1957) theory of superconduc- 
tivity Bardeen et al. (1959) have calculated the 
-lectronic contribution to the thermal conductivity 
when the dominant scatterers are impurities. Their 
-esult is similar to that given by previous theore- 
ical treatments by Heisenberg and Koppe 
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Abstract 


The thermal conductivity of niobium has been measured in the range 1.6 °K to 15 °K. Around 
the transition temperature (T,), thé main part of the heat was carried by the electrons, and the 
observed thermal conductivity in the Superconducting phase was found to be in good agreement 


Below T = 2.5 *K the thermal conductivity was dominated by the lattice contribution, and the 
data gave information about the dislocation density. A small lattice contribution to the total thermal 


(1950) and in fair agreement with experiments on 
Sn, Zn, Al and In. For a detailed account see the 
review by BardeenandSchrieffer (1961). 
The main purpose of the present paper is to com- 
pare some measurements of the thermal conduc- 
tivity of superconducting niobium with the above 
mentioned theory of Bardeen and co-workers. 

We have measured the thermal conductivity of 
a niobium sample in the range kK: GEEKY tombe "Ej 
and in the lower part of this range, below 2.5 °K, 
the thermal conductivity is dominated by the lattice 
part. Since the contributing phonons will then 
mainly be scattered by lattice defects, we have used 
the measured thermal conductivity to get informa- 
tion about the dislocation density. Furthermore, 
we have also observed a small lattice contribution 
to the thermal conductivity above Te. 

Following well established rules we will denote 
the lattice and electronic part of the thermal con- 
ductivity by the subscripts g and e respectively. 
When we refer to the superconducting phase we 
will use the subscript s and similarly » when we 
mean the normal phase. Hence we have for the 
total thermal conductivity 


K, = Ke, ER or K, = Kgs KES (1) 
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2. Sample and experimental procedure. 


The sample was supplied by Johnson, Matthey 
and Co. Ltd. and was in the form of a rod 10 cm 
long and 0.474 cm in diameter. The sample was 
measured as received without any annealing, and 
a spectroscopical test gave the impurities listed in 
table I. 


Table I. The impurities in the niobium sample deter- 
mined by a spectrographical test. 


Impurity Parts per million 
Ta About 10 0 
Fe 300 
Si 200 
Sn 50 
Pb 30 
Cu 5 
Ni 5 
Mg 1 


Fig. 1 shows the experimental cryostat schema- 
tically. The thermal conductivity was measured in 
a conventional way. A uniform temperature gra- 
dient was established in the rod (D) using a 
heating element (F), and the temperature was 
measured at two points along the rod (E). The 
cold end of the sample was attached to the heat 
reservoir (B). 

For measurements at temperatures below 4.2 °K 
the reservoir was filled with liquid helium and by 
pumping through a manostat system the tempe- 
rature could be kept constant within 0.001 deg at 
any stage down to 1.5 ^K. For measurements at 
temperatures above 4.2 *K a suitable temperature 
gradient was established along the tubes connec- 
ting the reservoir and the helium bath. This was 
obtained by a heating element wrapped around 
the reservoir (A). 

Two carbon resistors were used as thermome- 
ters. They were calibrated against the helium vapor 
pressure below 4.2 *K and against a platinum 
resistance thermometer (C) above 10 °K. The 
latter was calibrated in the range 107 to 90 *K by 
the National Bureau of Standards. We used the 
following interpolation curve (Clement and 
Quinell, 1952): 


A, + A,/T = log R + Azllog R (2) 


where the A’s are constants and R is the measured 
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Fig. 1. Experimental chamber, schematic. A-heating ele- 
ment, B-heat reservoir, C-calibrated resistance thermo- 
meter, D-sample, E-carbon resistor thermometers, and 
E-heating element. 1 


resistance. Several calibration points were taken 
and after determining the constants Ao, A, and 
A, the remaining points were used to plot devi- 
ations from (2) versus T. We estimate the uncer- 
tainty in the temperature determination to be less 
than 0.01 deg below 12 *K and somewhat larger 
above this temperature. i “i 


The thermal conductivity was determined from 
the relation 


Ke (3) 


here a is the cross section, Q the heat input and 
(A T the temperature difference between the two 
{measuring points separated by the distance /. We 
iestimate the uncertainty in K to be 5 per cent. 


3. Results and discussion. 

The measured thermal conductivity as function of 
the temperature is shown in fig. 2. The figure also 
shows some results from other investigators on 


T 7 1 
K E 


T 


[= 1 
۱0 15°K 


er 


Fig. 2. The thermal conductivity in niobium versus T. 
. The triangles are the average of three measurements. 
The curve marked MR represent the measurements of 
Mendelssohn and Rosenberg (1952) on a 
| polycrystalline sample. The curves marked CMR were 
| measured by Calverley et al. (1961) and Rowell 
| (1961) on single crystals subjected to various amounts 
| of cold work. 


| niobium. The measurements of Mendelssohn 

| and Rosenberg (1952) were done on a poly- 

١ crystalline sample and those of Calverley et al. 

| (1961) and Rowell (1961) on single crystals 
subjected to various amounts of cold work. 

` Assuming that the thermal conductivity will be 
dominated by the electronic part (Ks ~ Kes i} 


| 
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except when T/T, «1, we can compare our mea- 
surements with the theory of Bardeen ef al. 
(1959) by plotting the observed K,/K, vs. Tp 
We have done this in fig. 3. The thermal conduc- 
tivity in the normal state was determined using 
the relation K, = Lo T, where L is the Lorentz 
ratio and 1/a = 1.81 p ohm cm is the measured 
electrical resistance above Te. Due to experimental 
difficulties during the measurements of the elec- 
trical resistance, we could not determine the tran- 
sition temperature of the sample better than 
و‎ (90208 ne KC 

Le v y et al. (to be published) have done ultra- 
sonic attenuation studies in a niobium sample of a 
purity comparable to that of our sample. They 
found T, = 8.9 °K and determined the supercon- 
ducting energy gap to be 2 e, = 3.5 &T,. Hence 
we have drawn the theoretical curve in fig. 5 using 
this value for the energy gap, and we have set 
1 9.0 °K. Fig. 3 shows that the observed 


= SET TUI T T am T T- 

10 ٢ ^ 

0,9 kn 4 
o8 H 4 
GZ iS | 
06 |- 4 
TE 4 
04 F E 
03 f : 
02 = 
0,1 zi 


ټغ 1 
1,0 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0 


Fig. 3. The observed Ks/Ka versus T/T. compared with 
the theory of Bardeen eż al. (1959). The transition 
temperature T. — 9.0 °K, and the superconducting 
energy gap 2 zo = 3.5 Te. 


K,[K, has qualitatively the same temperature 
dependence as given by the theory. Since K, is the 
total thermal conductivity the deviation from the 
theory for small T/T, is obviously explained by 
the lattice contribution becoming dominant. The 
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somewhat slower fall of the measured value as 
compared to the theoretical value of K,/K, just 
below T, probably stems from the small lattice 
thermal conductivity included in the measured 
K,. As we will see later K, can be estimated to be 
about 7 per cent around T'e. 

From what is said above we find it safe to 
conclude that the theory of Bardeen e al. 
(1959) describe satisfactorily the electronic con- 
tribution to the thermal conductivity in niobium 
when the dominant scatterers are impurities. 

It is to be expected that the thermal conductivity 
in the lower part of the temperature range will be 
dominated by the lattice contribution, and given 
by an expression of the form 

1/K,, — W,,— A[T3 + B|T2 (4) 
Here the first term expresses the lattice heat resi- 
stivity (Ws) as limited by scattering with crystal 
boundaries and the latter term scattering due to 
dislocations. 

The measured thermal conductivity below 
2.5 *K seems to yield a power in T of approxi- 
mately 1.8 rather than a value larger than 2 as 


400- W ‘cm deg. 
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follows from (4). Although our observations are 
too few to say anything definite on the T-depen- 
dence it is interesting to note that Rowell 
(1961) finds a decreasing power in T with incre- 
asing dislocation density. He finds a T1-9-depen- 
dence for his most deformed specimen. This could 
possibly mean that another term in (4) is neces- 
sary to give a complete description of the lattice 
thermal conductivity in samples with large dislo- - 
cation densities. To estimate the dislocation den- 
sity in our sample we have fitted a curve of the - 
form K,,— T?/B to the observations in the lowest 
part of the temperature range as shown in fig. 4. - 
Following Klemens (1958) we have used the 
outcoming B to calculate the dislocation density 
N, and we find N = 5.2۰1012 cm. However, 
Lomer and Rosenberg (1959) have com- 
pared dislocation densities determined this way - 
with the densities found by electron microscopy, 
and find that Klemens’ formula yields dislo- - 
cation densities about an order of magnitude too 
large. 


Di 
20-10? Wem deg. 
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small lattice. contribution to the measured 
ermal conductivity can be observed above Te. 
nis is best illustrated in a linear plot of K versus 
as in fig. 5. Here we have drawn the line 
-» = Lo T as determined from the measured 
ectrical conductivity. We see from fig. 5 that all 
served values for K, fall above this line. Con- 
ering the rather impure sample in question, this 
most likely explained by the lattice term be- 
ming noticeable. We estimate K,, to be about 7 
t cent of K, around Te 
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